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Abstract: In this article we unravel the role of matter effect in neutrino oscillation in 
the presence of lepton-flavor-conserving, non-universal non-standard interactions (NSI’s) of 
the neutrino. Employing the Jacobi method, we derive approximate analytical expressions 
for the effective mass-squared differences and mixing angles in matter. It is shown that, 
within the effective mixing matrix, the Standard Model (SM) VE-exchange interaction only 
affects 012 and 0i3, while the flavor-diagonal NSI’s only affect 023- The CP-violating phase 
5 remains unaffected. Using our simple and compact analytical approximation, we study 
the impact of the flavor-diagonal NSI’s on the neutrino oscillation probabilities for various 
appearance and disappearance channels. At higher energies and longer baselines, it is found 
that the impact of the NSI’s can be significant in the —)• channel, which can probed 
in future atmospheric neutrino experiments, if the NSPs are of the order of their current 
upper bounds. Our analysis also enables us to explore the possible degeneracy between the 
octant of 023 and the sign of the NSI parameter for a given choice of mass hierarchy in a 
simple manner. 
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1 Introduction and Motivation 


The recent measurement of the moderately large value of the 1-3 mixing angle [1-11], 
quite close to its previous upper limit [12, 13], strongly validates the standard three-flavor 
oscillation model of neutrinos [14, 15], which has been quite successful in explaining all 
the neutrino oscillation data available so far [16-18], except for a few anomalies observed 
at very-short-baseline experiments [19]. This fairly large value of @13 greatly enhances the 
role of matter effects^ [21-23] in currently running and upcoming long-baseline [24-26] and 
atmospheric [27-30] neutrino oscillation experiments aimed at determining the remaining 
fundamental unknowns, in particular, the neutrino mass hierarchy,^ possible presence of 
a CP-violating phase <5, and the octant ambiguity of 623 [31] if the 2-3 mixing angle is 
non-maximal. A clear understanding of the sub-leading three-flavor effects [32, 33] in the 
neutrino oscillation probabilities in matter is mandatory to achieve the above goals. 

Furthermore, in addition to the Standard Model (SM) IT-exchange interaction, various 
models of physics beyond the SM predict non-standard interactions (NSI) of the neutrino 
[21, 34-37], which could affect neutrino propagation through matter. Indeed, such NSI’s^ 
arise naturally in many neutrino-mass models [40-55] which attempt to explain the small 
neutrino masses and the relatively large neutrino mixing angles, as suggested by current 
oscillation data [16-18], as well as in many other models to be discussed in a later section. 
Thus, understanding how the presence of NSI’s would affect the three-flavor neutrino os¬ 
cillation probabilities in matter is crucial in extracting the fundamental unknowns listed 
above, and also in searching for new-physics signatures in neutrino oscillation data. 

While the three-flavor oscillation probabilities in matter can be calculated numerically 
on a computer, with or without NSI’s, and properly taking into account the changing 
mass-density along the baseline, the computer program is a blackbox which does not offer 
any deep understanding as to why the probabilities depend on the input parameters in 
a particular way. If the mass-density along the baseline is approximated by an average 
constant value, then exact analytical expressions for the three-flavor oscillation probabilities 
can, in principle, be derived as was done for the SM case in Refs. [56-59]. But even for 
the SM case, the expressions are extremely lengthy and too complicated to yield much 
physical insight. Thus, for the SM case, further approximations which simplify the analytic 
expressions while maintaining the essential physics have been developed by various authors 
[60-71] to help us in this direction. Indeed, approximate analytical expressions of the SM 
neutrino oscillation probabilities in constant-density matter have played important roles in 
understanding the nature of the flavor transitions as functions of baseline L and/or neutrino 
energy E [65-67, 69]. To obtain similar insights for the NSI case, approximate analytical 
expressions for the three-flavor oscillation probabilities in constant-density matter in the 
presence of NSI’s are called for. Once they have provided us with the intuition we seek, 
on how and why the oscillation probabilities behave in a particular way, we can then 
resort to numerical techniques to further refine the analysis, e.g. taking into account the 

^For a recent review, see Ref. [20]. 

^There are two possibilities: it can be either ‘normal’ if = ml — m\ > 0, or ‘inverted’ if 5m|i < 0. 

^Present status and future prospects of NSI’s are discussed in recent reviews [38, 39]. 
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non-constant mass-density, if the need arises. 

In previous works [ 68 , 70], we looked at the matter effect on neutrino oscillation due to 
the SM Vk-exchange interaction between the matter electrons and the propagating electron 
neutrinos. Employing the Jacobi method [72], we showed that the said matter effects could 
be absorbed into the ‘running’ of the effective mass-squared-differences, and the effective 
mixing angles 612 and ^13 in matter as functions of the parameter a = 2 \/ 2 GFNeE, while 
the effective values of O23 and the CP-violating phase 6 remained unaffected. Here, Gp 
is the Fermi muon decay constant, Ne is the electron density, and E is the energy of 
the neutrino. The approximate neutrino oscillation probabilities were obtained by simply 
replacing the oscillation parameters in the vacuum expressions for the probabilities with 
their running in-matter counterparts. 

This running-effective-parameter approach has several advantages over other approaches 
which approximate the neutrino oscillation probabilities directly. First, the resulting ex¬ 
pressions for the probabilities are strictly positive, which is not always the case when the 
probabilities are directly expanded in some small parameter, and the series truncated af¬ 
ter a few terms. Second, the behavior of the oscillation probabilities as functions of L 
and E can be understood easily as due to the running of the oscillation parameters with 
a, as was shown in several examples in Refs. [ 68 , 70]. Third, as will be shown later, it is 
very convenient in exploring possible correlations and degeneracies among the mass-mixing 
parameters that may appear in matter in a non-trivial fashion. 

In this paper, we extend our previous analysis and investigate how our conclusions are 
modified in the presence of neutrino NSI’s of the form 

^NC-NSi = - 2\/2G'Fg^g [fliiPcf) , (1-1) 


where subscripts a, (3 = e, /r, r label the neutrino flavor, / = e,u,d mark the matter 

fC 

fermions, G = L,R denotes the chirality of the // current, and are dimensionless 
quantities which parametrize the strengths of the interactions. The hermiticity of the 
interaction demands 

dS = (d?)' ■ (1-2) 

For neutrino propagation through matter, the relevant combinations are 




- E 

f=e^u,d 


j ^ 


E 

f=e,u,d 


Ne ’ 


(1.3) 


where Nf denotes the density of fermion /. In this current work, we limit our investigation 
to flavor-diagonal NSI’s, that is, we only allow the with a = /3 to be non-zero. The 
case of flavor non-diagonal NSI’s will be considered in a separate work [73]. 

In the following, we will show how the presence of such flavor-diagonal NSI’s affect 
the running of the effective neutrino oscillation parameters (the mass-squared differences, 
mixing angles, and CP-violating phase), and ultimately how they alter the oscillation 
probabilities. We find that due to the expected smallness of the Saa's as compared to the 
SM IT-exchange interaction, there is a clear separation in the ranges of a at which the 
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NSI’s and the SM interaction are respectively relevant. Furthermore, within the effective 
neutrino mixing matrix, the SM interaction only affects the running of 612 and ^ 13 , while 
the flavor-diagonal NSI’s only affect the running of 023- The CP-violating phase 5 remains 
unaffected and maintains its vacuum value. 

We note that similar studies have been performed in the past by many authors e.g. 
in Refs. [69, 74-81]. This work differs from these existing works in the use of the Jacobi 
method [72] to derive compact analytical formulae for the running effective mass-squared 
differences and effective mixing angles, which provide a clear and simple picture of how 
neutrino NSI’s affect neutrino oscillation. We also note that we have addressed the same 
problem with a similar approach previously in Refs. [82-85]. The current paper updates 
these works by allowing for a non-maximal value of 622 ,, a generic value of the CP-violating 
phase (5, a larger range of a, and refinements on how the matter effect is absorbed into the 
running parameters. 

This paper is organized as follows. We start section 2 with a discussion on neutrino 
NSI’s: how and where they arise, how they affect the propagation of the neutrinos in matter, 
and show that the linear combinations relevant for neutrino oscillation are rj = err)/2 

and C, = £ee — + £Tr)/2. This is followed by a brief discussion on the theoretical 

expectation for the sizes of these parameters, and their current experimental bounds. In 
section 3, we use the Jacobi method to calculate how the NSI parameter rj affects the 
running of the effective mass-squared differences, effective mixing angles, and the effective 
CP-violating phase as functions of a = a(l-|-C) for the neutrinos. To check the accuracy of 
our method, we also present a comparison between our approximate analytical probability 
expressions and exact numerical calculations (for constant matter density) towards the end 
of this section. Section 4 describes the advantages of our analytical probability expressions 
to figure out the suitable testbeds to probe these NSPs of the neutrino. In this section, we 
also present simple and compact analytical expressions exposing the possible correlations 
and degeneracies between 023 and the NSI parameter r/ under such situations. Finally, we 
summarize and draw our conclusions in section 5. The derivation of the running oscillation 
parameters for the anti-neutrino case is relegated to appendix A where we also compare 
our analytical results with exact numerical probabilities. In appendix B, we examine the 
differences in the exact numerical probabilities with line-averaged constant Earth density 
and varying Earth density profile for 8770 km and 10000 km baselines. 

2 A Brief Tour of Non-Standard Interactions of the Neutrino 

In this section, we first briefly discuss the various categories of models that give rise to 
NSFs of neutrino. 

2.1 Models that Predict NSI’s of the Neutrino 

NSFs of the neutrino arise in a variety of beyond the Standard Model (BSM) scenarios 
[38, 82-86] via the direct tree-level exchange of new particles, or via flavor distinguishing 
radiative corrections to the Zvv vertex [87-89], or indirectly via the non-unitarity of the 
lepton mixing matrix [90]. BSM models which predict neutrino NSFs include: 
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1. Models with a generation distinguishing Z' boson. This class includes gauged L^ — L^ 
and gauged Lg — L-r [91, 92], gauged B — aL^ — /3L^ — ^Lr (with a + /3 + 7 = 3) 
[89, 93-97], and topcolor assisted technicolor [87, 98]. 

2. Models with leptoquarks [99-103] and/or bileptons [104]. These can be either scalar 
or vector particles. This class includes various Grand Unification Theory (GUT) 
models and extended technicolor (ETC) [105-107]. 

3. The Supersymmetric Standard Model with ii-parity violation [ 88 , 108-111]. The 
super-partners of the SM particles play the role of the leptoquarks and bileptons of 
class 2 . 

4. Extended Higgs models. This class includes the Zee model [52], the Zee-Babu model 
[53, 112, 113], and various models with SU{2) singlet [ 86 , 114-117] and SU{2) triplet 
Higgses [115, 118, 119], as well as the generation distinguishing Z' models listed under 
class 1 . 

5. Models with non-unitary neutrino mixing matrices [19, 90, 120-134]. Apparent non- 
unitarity of the mixing matrix for the three light neutrino flavors would result when 
it is part of a larger mixing matrix involving heavier and/or sterile fields. 

Systematic studies of how NSI’s can arise in these, and other BSM theories can be found, for 
instance, in Refs. [38, 85, 86 , 135]. Thus, the ability to detect NSTs in neutrino experiments 
would complement the direct searches for new particles at the LHC for a variety of BSM 
models. Next, we focus our attention to see the role of lepton-flavor-conserving NSI’s when 
neutrinos travel through the matter. 


2.2 Lepton-Flavor-Conserving NSI’s 

The NSI’s of Eq. (1.3) modify the effective Hamiltonian for neutrino propagation in matter 


in the flavor basis to 



( 

9 

rui 

0 

0 ’ 


U 

0 

2 

mi 

0 

+ a 

2E 
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0 

0 

ml_ 



1 + £, 


ee ^e/i ^er 

err_ 


( 2 . 1 ) 


where U is the vacuum Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix [136-138], E 
is the neutrino energy, the matter-effect parameter a is given by 


= 2V2GFNeE = 7.6324 x 10"^(eV^) 


P 


g/cm^ 


A) 


For Earth matter, we can assume ~ Np = Ag, in which case 

_ _ I Q 

ta/3 — ^ ^0/3 + '3^0/3 


+ 3^0/3 


( 2 . 2 ) 

Nd K, 3Ne. Therefore, 
(2.3) 


Since we restrict our attention to lepton-flavor-conserving NSI’s in this paper, the 
effective Hamiltonian in the flavor basis takes the form 
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mf 0 0 

0 ml 0 

+ a 

V 
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1 + £ee 0 I 

0 Epfj, I 
0 0 


(2.4) 
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In the absence of off-diagonal terms, we can rewrite the matter effect matrix as follows: 


1 + £ee 0 0 

0 0 

0 0 Err 

1 , ^ 

J- “T c-ee 2 

= 0 

0 


0 0 




= o 





+ Ej 


1 0 

0 0 


-|- O(e^) terms -|- unit matrix term , 



1 0 0 
0 1 0 
0 0 1 


(2.5) 


where we have assumed that the e’s are small compared to 1. Let 


^TT 


C — ^ee 


“1“ ^TT 


( 2 . 6 ) 


and 

a = a(l + C). (2.7) 

Note that the NSI parameter r] introduced here is related to the parameter ^ which was 
used in Ref. [82] by 


^ = - 2 - 

Thus, the effective Hamiltonian can be taken to be 
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( 2 . 8 ) 


(2.9) 


Note that the rj = 0 case simply replaces a with a = a(l -|- () in the SM Hamiltonian. 
Thus, the effective mass-squared differences, mixing angles, and CP-violating phase have 
the same functional dependence on a as they had on a in the SM case. In other words, a 
non-zero C simply rescales the value of a by a constant factor. The presence of a non-zero 
r], on the other hand, requires us to diagonalize H with a different unitary matrix from 
the r] = 0 case, and this will introduce corrections to the effective oscillation parameters 
beyond a simple rescaling of a. Next, we discuss the constraints that we have at present 
on these NSI parameters. 


2.3 Existing Bounds on the NSI parameters 

Before looking at the matter effect due to the neutrino NSI’s, let us first look at what is 
currently known about the sizes of the and the combinations rj = (e^^ — ErT)l‘^ and 

C — ^ee 
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2.3.1 Theoretical Expectation 

Theoretically, the sizes of the NSI’s are generically expected to be small since they are 

putatively due to BSM physics at a much higher scale than the electroweak scale, or to 

loop effects. If they arise from the tree level exchange of new particles of mass A, which 

fC 

would be described by dimension six operators, we can expect the ’s to be of order 
0(Af^/A^). Loop effects involving a heavy particle of mass A would be further suppressed 
by a factor of 0(l/47r) or more. Processes that lead to dimension eight operators would 
be of order 0(M^/A‘^). Thus, if we assume A = 0(1 TeV), we expect the e^^’s, and 
consequently rj and C; to be 0(10“^) or smaller. 

2.3.2 Direct Experimental Bonnds 

Direct experimental bounds on the flavor-diagonal NSI parameters ’s are available from 
a variety of sources. These include t'eC scattering data from LAMPF [139] and LSND 
[140], scattering data from the reactor experiments Irvine [141], Krasnoyarsk [142], 
Rovno [143], MUNU [144], and Texono [145], Uf,q scattering data from CHARM [146], 
scattering data from CHARM H [147, 148], scattering data from NuTeV [149- 
151], e"’'e“ —)• zzpy data from the LEP experiments ALEPH [152-154], L3 [155-157] OPAL 
[158-161], and DELPHI [162], and neutrino oscillation data from Super-Kamiokande [163], 
IceCube and DeepCore [164, 165], KamLAND [166], SNO [167], and Borexino [168, 169]. 
These have been analyzed by various authors in Refs. [170-192], and collecting the results 
of the most recent analyses, we place the following 90% C.L. bounds on the flavor-diagonal 
vectorial NSI couplings: 

kLI < 0.1 , le“el < 1 , Jefel < 1 , 

< 0.04 , < 0.04 , < 0.04 , (2.10) 

< 0.6 , 14^1 < 0.05 , le^.,1 < 0.05 . 

Note that the current experimental bounds on the NSI couplings are weak compared 
to the theoretical expectation of 0(10“^), the 90% C.L. upper bound on their absolute 
values ranging from 0(10“^) to 0(1). To combine these bounds into bounds on the Safs’s, 
we follow the procedure of Ref. [183], 

|e.?l S + , (2.11) 

and find 

\eee\ < 4 , < 0.2 , < 0.6 . (2.12) 

Neglecting possible correlations among these parameters, these bounds can again be com¬ 
bined to yield 

\v\ < \J = 0.3, 

ICI < ^keek + ~ ^ ' (2-13) 
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A tighter bound exists for 77 which has been obtained directly using solar and atmo¬ 
spheric neutrino data in Ref. [174], and more recently in Ref. [193] using atmospheric and 
MINOS data. In Ref. [174], only NSI’s with the d-quarks were considered and the following 
3 it (99.7% C.L.) bounds were obtained^: 

- 0.03 < < 0.02 , 

\et-eU < 0.05. (2.14) 

Since Nd = Nu = 3Ae, Ref. [174] is actually constraining Eapl^-, so this result can be 
interpreted as 


— 0.09 < e^r < 0.06 , 
|e^ 

Rescaling to 1.64cr (90% C.L.), we hnd 


^tt\ ^ 0.15 . 


Id! = 


^TT 


< 0.04 


Ref. [193] gives slightly different 90% C.L. bounds of 

< 0.035 , 

£tt| ^ 0.11 , 


(2.15) 


(2.16) 


(2.17) 


which translates to 



< 0.055 . 


(2.18) 


Thus, though r] and C are expected theoretically to be 0(10 ^), their current 90% C.L. 
experimental bounds are respectively ~0.05 and 0(1). 


3 Effective Mixing Angles and Effective Mass-Squared Differences 
— Neutrino Case 


3.1 Setup of the Problem 

As we have seen, in the presence of non-zero r] and C, the effective Hamiltonian (times 2E) 
for neutrino propagation in Earth matter is given by 


= U 

Ai 0 0 ’ 
0 A 2 0 

U = U 

’00 0 ’ 

0 Sm2i 0 

+ a 

'1 0 0 ’ 
0 0 0 

-bar/ 

’ 000 ’ 
0 10 


0 0 As 


0 0 


0 0 0 

✓ s. 

1 

0 

1 

1 —‘ 

1 _ 


= Ma ^ 

= Ho 


where a = a(l-l-C). The problem is to diagonalize Hr^ = Ho+arjMn and find the eigenvalues 
Aj {i = 1, 2, 3) and the diagonalization matrix U as functions of a and r], 

^The notation used in Ref. [174] is e = ej),- and e' = £% — 
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Parameter 

Best-fit Value & lu Range 

Benchmark Value 

dm|i 

(7.50 ±0.185) X 10-5 eV^ 

7.50 X 10-5 eV^ 

dmii 

(2.47l[];|]^?) X 10-5 eV^ 

2.47 X 10-5 eV^ 

sin^ 6*23 
023/° 

n 41+0.037 ^ Q rq-l-0.021 
’-’•^-‘■-0.025 T U.oy_o.022 

40.0l?;^ © 50.4+}J 

0.41 

023/rad 

0.698«©0.880tn2i 


sin^ 012 

0.30 ±0.013 

0.30 

012/° 

33.3 ±0.8 


012/rad 

0.580 ±0.014 


sin^ 013 

0.023 ±0.0023 

0.023 

013/° 

Q ^+0.44 
^•*^-0.46 


013 /rad 

0.15 ± 0.01 


5 /° 

J/tt 

3001^8 

1.67;0;57 

0 


Table 1. Second column shows the best-fit values and la uncertainties on the oscillation parameters 
taken from Ref. [194]. We use the values listed in the third column as benchmark values for which 
we calculate our oscillation probabilities in this work. 


To this end, we utilize the method used in Refs. [ 68 , 70] where approximate expressions 
for the Aj’s and U were derived for Hq, the rj = 0 case, using the Jacobi method [72]. The 
Jacobi method entails diagonalizing 2x2 submatrices of a matrix in the order which 
requires the the largest rotation angles until the off-diagonal elements are negligibly small. 
In the case of Hq, it was discovered in Refs. [ 68 , 70] that two 2x2 rotations were sufficient 
to render it approximately diagonal, and that these two rotation angles could be absorbed 
into ‘running’ values of 612 and @ 13 . The procedure that we use in the following for is 
to add on the term to Hq after it is approximately diagonalized, and then proceed 

with a third 2 x 2 rotation to rotate away the additional off-diagonal terms. 

As the order parameter to evaluate the sizes of the off-diagonal elements, we use 




V 

and consider Hq and H^^ to be approximately diagonalized when the rotation angles required 
for further diagonalization are of order = 0.005 or smaller. Note that we are using a 
slightly different epsilon (e) here to distinguish this quantity from the NSTs (£ 0 / 3 ). 

The eigenvalues A* (i = 1, 2, 3) and the diagonalization matrix U of H^ are necessarily 
functions of d = 2y/2GFNf.E{l + Q). In order to parametrize the size of d, we find it 
convenient to introduce the log-scale variable [ 68 , 70] 


so that d = corresponds to /3 = —2, and d = corresponds to /3 = 0. In the 

following, various quantities will be plotted as functions of (3. 
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Unless otherwise stated, we use the benchmark values of the various oscillation param¬ 
eters as given in the third column of Table 1 to draw our plots. These values are taken 
from Ref. [194] and correspond to the case in which reactor fluxes have been left free in the 
fit and short-baseline reactor data with L < 100 m are included. For sin^ 023) we consider 
the benchmark value which lies in the lower octant and CP-violating phase 5 is assumed to 
be zero. These choices of the oscillation parameters are well within their 3cr allowed ranges 
which are obtained in recent global ht analyses [16-18]. We also present results considering 
other allowed values of sin^ 023 and 6 which we discuss in section 3.5. In the evaluation of 
the sizes of the elements of the effective Hamiltonian, we will assume 0i3 ~ 0.15 = 0(e), 
cos(20i 2)/2 Ri 0.2 = 0(e), and |cos(2023)| ~ 0.18 = 0(e). We also assume that the NSI 
parameter ry is of order e^ = 0.03 (or smaller), since the current 90% C.L. (1.64cj) upper 
bound on |r/| was ~ 0.05, cf. Eqs. (2.16) and (2.18), though we allow it to be as large as 
0.1 in our plots to enhance and make visible the effect of a non-zero rj. 


3.2 Diagonalization of the Effective Hamiltonian 
3.2.1 Change to the Mass Eigenbasis in Vacuum 


Introducing the matrix 

Q 3 = diag(l,l,e*'^) , (3.4) 


we begin by partially diagonalizing the Hamiltonian as 


H' = qIu^HMQ3 


’00 0 ’ 


’100’ 


’000’ 

0 

0 

+ aQl[/'^ 

0 0 0 

ug3 +a0 

0 1 0 

0 0 


0 0 0 


0 0-1 


UQs, (3.5) 


Mn 


Mr, 


= M'aiOu, 013, 023) 


= M'(012, 013, 023, <5) 




where 


^a(012,013,023) — <33 

U^iUel U:,Ue3 


r c 2 

''12%3 

Cl 2 'Sl 2 cf 3 C 12 C 13 S 13 

UhUel Ul^Ue2 
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Cl 2 'Sl 2 cf 3 

'*12^13 'S12C13'S13 


U:,Uel U:,Ue2 U:,Ue3 


C 12 C 13 SI 3 

Sl2Cl3'Sl3 sfg 


(3.6) 

and 


M'(012, 013, 023, <5) 


<31 
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u;,u^2 u;,u,3 


'U;^Url U;^Ur2 U:^Ur3' 



u;2U,i u;,u^2 u;,u^3 

- 

U;^Url U:^Ur2 U*^Ur3 


V 

u;3Uf.i u;,u,2 u;^u,3 


u;^Uri u;,Ur2 u;,Ur3 
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sin( 20 i 2 ) sin(2023)si3 cos S + (^12 - C?2Si3) COs(2023) 




- Cj 2 e“')si 3 sin(2023) - (1 + Si 3 )s 12 C 12 cos(2023) 
-S 12 C 13 sin(2023)e“*^ C 12 S 13 C 13 cos(2023) 


> <33 


10 
























- cf 2 e *'^)si3 sin(26'23) - (1 + S? 3 )si 2 Cl 2 COs(26'23) 
- sin(20i2) sin(26*23)si3 cos 6 + ( 0^2 - ■s? 2 'Si 3 ) cos(2023) 
C 12 C 13 sin(26»23)e“*'^ + S 12 S 13 C 13 008(2023) 

-S12C13 sin(2023)e+*'^ + C12S13C13 cos(2023) 
C12C13 8111(2023)6+**^ + S12S13C13 COs(2023) 

-cfg COs(2023) 


(3.7) 


Using cos( 2023) = 0(e) and 0i3 = 0(e), we estimate the sizes of the elements of Ma to be: 


Ma 


0(1) 0(1) 0(e) 
0(1) 0(1) 0(e) , 

0(e) 0(e) 0(e2) 


(3.8) 


and those of Mr^ to be 




0(e) 0(e) 0(1) 
0(e) 0(e) 0(1) 
0(1) 0(1) 0(e) 


(3.9) 


Given that we have assumed r] = O(e^) or smaller, the off-diagonal elements of 
are suppressed compared to those of Ma, and only become important for a > \5m\i\, or 
equivalently, /? > 0. 


3.2.2 Diagonalization of a 2 x 2 hermitian matrix 

The Jacobi method entails diagonalizing 2x2 submatrices repeatedly. For this, it is 
convenient to note that given a 2 x 2 hermitian matrix. 


A 5e*^' 


1 0 


A B 


1 0 

B e-*^ C 


9 

T 

o 

_1 


B C 


0 e*^ 


the unitary matrix which diagonalizes this is given by 


2l, 5,0,0 e M , 


where 


That is. 


where 


Cuj 


1 0 


C-u) 


1- 

O 

T—1 

1_ 

-s^e“*+’ Coj 


Q 

T 

o 

_1 


C-ijj 


0 6*+* 


Co; = coso; , 


So) = sinw , 


tan 2a; 


2B 

C-A ■ 


t/t 


A 

Be-^^ C 


U 


Ai 0 
0 A2 


Acl - Csl ^ {A + C)t V{A - Cf + 4B^ 

4 - 4 2_ 

Oc^ - ^ {A + C)± V(A - 0)2 + 4B^ 

„2 _ 9 


(3.10) 


(3.11) 


(3.12) 


(3.13) 


(3.14) 
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(a) (b) 

Figure 1. (a) The dependence of on /3 = — log^ {a/\5m\i'\). (b) The /3-dependence of Aj_. 



(a) (b) 

Figure 2. (a) The dependence of s '^2 = sin0^2 ^od c ^2 = 0080^2 on /3 = — log^ {a/\5m\^\Y (b) 
The dependence of as \2 and ac'-y 2 on /3. The values are given in units of |i5to 3]^|. The asymptotic 
value of aci 2 is dm2iSi2Ci2/ci3 « 0.014 |(5m§i| = 0(e^|(5m§^|). 


the double signs corresponding to the two possible quadrants for 2a; that satisfy Eq. (3.12). 
In applying the above formula to our problem, care needs to be taken to choose the correct 
quadrant and sign combination so that the resulting effective mixing angles and mass- 
squared eigenvalues run smoothly from their vacuum values. 

3.2.3 r] = 0 Case, First and Second Rotations 

As mentioned above, we will first approximately diagonalize Hq, and add on the ar]M'^ term 
afterwards. Here, we reproduce how the Jacobi method was applied to in Refs. [68, 70]. 
There, a (1, 2) rotation was applied to Hq, followed by a (2,3) rotation, which was sufficient 
to approximately diagonalize Hq. 

1. First Rotation 
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Define the matrix V as: 


V = 

where 

Cip = cos (p , = sin if , tan 2 ip = 

Then, 

H'' = V^H'^V = 


Sip 0 

sp Cp 0 

0 0 1 


(3.15) 


ac^3 sin 26i2 


(5m|i — ac^g cos 26i2 


0<^<f 


OC 32 C 13 S 13 
J 


A'_ 0 

0 A'_)_ as'i2Ci3Si3 

ac'i2Ci'iSiz as'i2Ci2,si'i Srn^i + OS 33 


where 


and 


c'i2 = cos e'12 , S12 = sin e'i2 , 0'i2 = Gi 2 + P> , 


A± = 


(acfg + 6ml^) ± V(ac^3 - + 4ac^3S^2'^"i2 


21 


The angle 0^2 = ^12 + V can be calculated directly without calculating ip via 

5m2i sin 2612 


tan 26 'i 2 = 


^12 <0'i2<-] ■ 


577123 cos 2012 — hcfg ’ 

As j3 is increased beyond —2, the A'_i_ asymptote to 

A+ —)• aci3 + 5m2iSi2 , 
\'_ —)• 57T72iCi2 • 


- ^12^ 

(3.16) 

(3.17) 

(3.18) 

(3.19) 

(3.20) 


(3.21) 


The dependences of 0(2 ^^cl A'_i_ on (5 are plotted in Fig. 1. Note that 0(2 increases 
monotonically from 0i2 to 7r/2 with increasing (5. The /3-dependence of s'12 = sin 0^2 
and C32 = cos 032 are shown in Fig. 2(a). As /3 is increased beyond —2, that is 
a = 5?7i23, S32 grows rapidly to one while C32 damps quickly to zero. In fact, the 
product ac32 stops increasing at around [3 = —2 and plateaus to the asymptotic 
value of 5m23Si2Ci2/c33 ss 0.014 |5?7733| = 0(e^)|5m33| as shown in Fig. 2(b). That 
is: 


OS 32 = |5?T7|3|0(e ^) , 

ac32 = |5?T7|3| 0(e“™™^^’“^^) < |5?77|3| O(e^) . (3.22) 

Note also that the scales of A'_j_ are given by 

X'j^ = 0(max(5?T723, a)) = , 

A'_ = 0(min(5m2i,a)) = |5?77|3| 0(e“. (3.23) 
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p 


3 n 
8 


0 


8 



P 


(a) 


(b) 


Figure 3. The dependence of (a) cj) and (b) = 6*13 + (/> on /3 = — log^ (a/|(5m§^|) for the normal 

(NH) and inverted (IH) mass hierarchies. 


2. Second Rotation 

Given that the (1, 3) element of Hq, namely dc'^ 2 Ci 3 'Si 3 , is at most of order \5m\i\0{e^) 
for all d, whereas the ( 2 ,3) element ds']^ 2 Ci 3 'Si 3 will continue to increase with d, it 
is the (2,3) submatrix that needs to be diagonalized next. The matrix W which 
diagonalizes the ( 2 ,3) submatrix of Hq is 


W = 


1 0 0 
0 
0 


(3.24) 


where = cos (/>, = sin cj), and 

, ^ 2 ds; 2 'Si 3 Ci 3 _ _ a sin 26>i3 _ 

(5m|i + dsfg — A'_)_ (<^"^-31 ~ ^'itT‘ 21 ^ 12 ) ~ d cos 2^13 

The angle 4> is in the first quadrant when > 0 (normal hierarchy) and increases 
from zero to ^ — 0 i 3 as /? is increased, (j) is in the fourth quadrant when < 0 

(inverted hierarchy) and decreases from zero to —^13 as /3 is increased. This in¬ 
dependence of (j) is shown in Fig. 3(a) for both mass hierarchies. 

Using W, we obtain 




W^H'qW 

A'_ —dC32Cl3Si3S0 dC32Cl3Si3C,^ 

dCi2Ci3Si3S(j() A^ 0 

dC32Cl3Si3C0 0 A'_t 


(3.26) 


where the upper signs are for the > 0 (normal hierarchy) case and the lower 

signs are for the < 0 (inverted hierarchy) case, with 


[A'+ + (dm|i -b dsfg)] ± y/[A'+ - {6ml^ + dsfg)]^ -b 4(ds'32Ci3Si3)^ 


(3.27) 
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(a) Normal Hierarchy 


(b) Inverted Hierarchy 


Figure 4. The /3-dependence of for the (a) normal and (b) inverted mass hierarchies. 


As j3 is increased beyond 0, the asymptote to 

A” —^ d -|- 6rn^iSi^ + ^^ 21 ^ 12^^13 > 

A” —)• -|- ^^21^12^13 ) (3.28) 

for both mass hierarchies. Note that A(( < 0 for the < 0 case. The in¬ 

dependences of Aj. are shown in Fig. 4. Order-of-magnitude-wise, we have 

A" = , A'l = , if > 0 , 

A" = , |A"| = |<5mii|0(l) , if < 0 . 

(3.29) 

In particular, in the range /3 > 0 we have 

\"_ = 0(1) , A" = \5m\i\0{e~^) 

A'l = |dm|i|0(e-^) , |A'^| = |dm|i|0(l), 

For the off-diagonal terms, since dc'^2 = l'^"^3i|0(e^), C13 
0(l/e), C0 = 0(e/l), we have 

- ac'i2Ci3Si3S,^ = 0(e^/e^), 
ac'i2Ci3Si3C,^ = 0(e^/e^). 


, if > 0 , 

if < 0 . 


(3.30) 


= 0(1), si 3 = 0(e), = 


(3.31) 


Thus, looking at the sizes of the elements of Hq in that range we find: 


tt/// 

^0 




0(e2) 0(e3/e^) O(eVe^) 

0(e3/e^) 0(l/e-^) 0 

0(eVe3) 0 0{e-^/l) 


(3.32) 


where the elements with two entries denote the two different mass hierarchies, 0(NH/IH), 
and we can see that further diagonalization only require angles of order O(e^). There¬ 
fore, Hq can be considered approximately diagonal. 
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3.2.4 r] ^ 0 Case, Third Rotation 

Let us now consider the 7 ^ 0 case. If we perform the same (1,2) rotation V on = 
Hq + ar]M'^ as we did on Hq, the M' part is transformed to 

M^{9i2,6i3,623,6)V = (012^-1^, 013, 023, 5) 

= ^'12 

= M;(0i2, 013, 023,-5) . (3.33) 


Using 0(2 —t f, asi 2 0{e^)\6m‘^^\ as j3 is increased beyond —2, we can 

approximate 


ar]M'^ (012,013,023,5) 

~ a7?M' (^|, 013, 023, 5 ) 


= at] 


cos( 2023) e*'^si3 sin(2023) 

e“*'^si3sin(2023) -sfg cos(2023) 
-e“*^ci3 sin(2023) S 13 C 13 cos(2023) 


-e*‘^ci3 sin(2023) 
S13C13 cos(2023) 
-cfg COs(2023) 


Performing the (2, 3) rotation W next, we find: 

iptM; (^^,013,023,5) vu 
= -^(7 ( ^ \013 + (^ , 023, 5 ) 

= 0(3 

= M; (^,013,023,5) 

cos( 2023) e*'^Si 3 sin(2023) -e^'^-^is sin(2023) 

= e"*'^S 33 Sin( 2023 ) -S13 cos(2023) s'lgc'ig cos(2023) 

-e“*'^C33sin(2023) 513013003 ( 2023 ) -c'^ cos(2023) 


(3.34) 


(3.35) 


where s '13 = sin 0 'i 3 and c '13 = cos 0 'i 3 . The angle 0 'i 3 = 0 i 3 + (p can be calculated directly 
without the need to calculate (p using 


tan 20(3 


((5m|i — (5m|isf2) sm20i3 
((5m|i — hm|iSi 2 ) cos 20 i 3 — a ’ 


(3.36) 


and its /3-dependence is shown in Fig. 3(b). As can be seen, 0 'i 3 increases rapidly to 7 r /2 
when 5m|i > 0 , while damping quickly to zero when 5m|i < 0 , once /3 is increased above 
zero. Consequently, acos 0 'i 3 for the 5m|i > 0 case, and asin 0 'i 3 for the 5m|i < 0 case 
plateau to 013513(1 — e^'Si 2 )| 5 m|i| = 0{e)\5m\i\ as /3 is increased as shown in Fig. 5. Note 
that in the 5m,|i > 0 case, a cos 0^3 increases to 0{l)\5m\i\ in the vicinity of /3 = 0 before 
plateauing to 0 (e)| 5 m 3 i|. This will cause a slight problem in our approximation later. We 
now look at the normal and inverted mass hierarchy cases separately. 


1. ^mgi > 0 Case 
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(a) Normal Hierarchy 



P 

(b) Inverted Hierarchy 




P 


P 


(c) Normal Hierarchy 


(d) Inverted Hierarchy 


Figure 5. The dependence of sindj^g and on /3 = — log^ {a/\5m\i\') for the (a) normal and 

(b) inverted mass hierarchies. The dependence of asin0gg and a cos 6*^3 on /3 = — log^ (a/l^rnggl) 
for the (c) normal and (d) inverted mass hierarchies. 


For the > 0 case dcgg —)• 0{e)\5m\i\ as /3 is increased beyond 0. Therefore, we 
can approximate 


H'” = WW^H'^VW 

« < + a7yM;(|, 013,023, 5 ) 

A'_ + d?7 cos( 2023) 0776**^513 sin(2023 ) 0 

~ dr7e“*^sl3 sin(2023) A" — dr/sl3 cos(2023 ) 0 
0 0 A'l 


(3.37) 


where we have dropped off-diagonal terms of order 0(e^)|(5777|g| or smaller. (This 
approximation breaks down in the vicinity of /3 = 0 where both dsl 3 and dcl 3 are of 
order 0{l)\6m\i\.) Define the matrix X as 


X = 


s^e''^ 0 

— 5^6“**^ 0 

0 0 1 


(3.38) 
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where c-y = cos x, ■Sy = sin x, and 


tan2x = 


2ar]s[^ sin(26*23) 


(A" - A'_) - ari{l + 003(26*23) 

2 ar] sin( 2023 ) 

- sf 2 S? 3 )] “ 2 a 7 y008(26*23) ' 


Note that 


vr 


0 < X < 2 “ ^23 for 7/ > 0 , - 6*23 < X < 0 for r/ < 0 . 


(3.39) 


(3.40) 


The / 3 -dependence of x is shown in Fig. 6 for several values of 77, both positive 
(Fig. 6(a)) and negative (Fig. 6(b)). 


Using X, we find 


tthh _ vt ij'" Y 
— -A. -^Tj 


\'x_ 0 0 

0 A'^, 0 
0 0 A'l 


(3.41) 


where 


X'" 


(A" + A'_ -|- 0770^3 cos 26*23) i ^ [A" — A'_ — 077(1 + S33) cos 2^23]^ + 4(077533 sin 26*23)^ 


(3.42) 


Thus, is approximately diagonal. The asymptotic forms of Xx± at /3 S> 0 are 


A 


/// 

x+ 


A 


/// 

X- 


a\r]\ + 
-a\v\ + 


<5"l3lC?3«23 + <^"i2l(c?2ci3 + «12S?3S23) 
+ '5"^2i(c?2S23 + 

+ '5"^il(c?2S23 + 

+ 5ml^{4243 + «12S?3«23) 


for 77 > 0 
for 77 < 0 

for 77 > 0 
for 77 < 0 


(3.43) 


The / 3 -dependence of Xx± are shown in Figs. 6(c) to 6(f). 


2. Srn^^ < 0 Case 

For the <577733 < 0 case we have as'^g —)• 0 (e)|<577733 1 as /3 is increased beyond 0 . 
Therefore, we can approximate 


H'” = W^V^H'^VW 


H”' + 


0 -r ^) ^13) ^23) <5 

A'_-I-077003(26*23) 0 —a77e*'^C33 sin( 2023 ) 

0 A" 0 

-a 77 e“*'ic 33 sin( 2023 ) 0 A'/. — 077033 cos( 2023 ) 


(3.44) 


where we have dropped off-diagonal terms of order 0(e^)|<577733| or smaller. Unlike 
the 577733 > 0 case, this approximation is valid in the vicinity of /3 = 0 since 0533 
never exceeds 0(e) 15777331 for all o. 
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(a) 7? > 0 


(b) ?7 < 0 



yS 

(e) 77 > 0 (f) r; < 0 


Figure 6 . The / 3 -dependence of x and X'xj. for several values of 77 with S23 = 0 . 4 . 


Define the matrix Y as 


Y 


C'lp 0 SijjG 
0 10 

Sipc 0 


(3.45) 
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where = cos ip, = simp, and 


tan 2^ = 


2ar](P^^ sin(2023) 

(A" - A'_) - a?7(l + c'^g) cos(2023) 

2 ar] sin(2023) 

[|(5m|i|c?3 + 5m|i(cf2 - 'Si2'Si3)] + 2ar] cos{2e23) ’ 


Note that 


— 6*23 < V’ < 0 for r/ < 0 , 


0 < 'tp < — — 023 for r/ > 0 . 


(3.46) 


(3.47) 


Comparing Eq. (3.39) and Eq. (3.46), we can infer that 'ip^'t]) ~ xi~'n)^ the small 
difference due to the Sm^i term in the denominator of the expressions for tan 2 x and 
tan2'i/’. This can be seen in Eig. 7 where the /^-dependence of ip is shown for several 
values of rj, both positive (Eig. 7(a)) and negative (Eig. 7(b)). 


Using Y, we find 


= Y^H"'Y 


where, 


A"'+ 0 0 

0 A" 0 

0 0 A'^'. 


(3.48) 


A 


/// 

Y± 

{X'!_ + X'_ + arjs'i^ cos 


26 * 23 ) ± ^[A" 


A'_ - ari{l + cd) cos 26 * 23 ]^ + 4 (arycg 3 sin 26 * 23 )^ 

2 ■ 

(3.49) 


Thus, H”P_ is approximately diagonal. The asymptotic forms of Ay^ at /? S> 0 are 


Ayy —)> a| 7 y| -I- 

Ay_ —> + 


-|^"i3l|Cl3«23 + <5"l2l(Cl2C23 + Sl2«13«23) 
-|^"i3l|Cl3C23 + '^"i2l(Cl2S23 + 512^13^23) 

-|(5m|i|cf3ci3 + dml^icl^sh + 

-|<5"l3l|c?3«23 + '5w2l(c?2ci3 + S?2Sl3«23) 


for T] > 0 
for r] < 0 

for r] > 0 
for r] < 0 


(3.50) 


The /3-dependence of Ayj_ are shown in Eig. 7(c) to Eig. 7(f). 


3.3 Effective Mixing Angles for Nentrinos 

We have discovered that the unitary matrix which approximately diagonalizes is U = 
UQ 3 VWX when > 0, and U = UQ 3 VWY when 5rn^^ < 0. Introducing the notation 



cos 6 

sin 6 e 

0 ’ 

Ri2{0,5) = 

— sin 6* 


cos 9 

0 


0 


0 

1 


cos 6 

0 

sin 9 e~ 

i5 

7 / 13 ( 0 , (3) = 

0 

1 

0 



— sin 6* 

0 

cos 9 
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(c) rj > 0 


(d) r? < 0 




yS 




(e) rj > 0 


(f) r? < 0 


Figure 7. The ^-dependence of ip and Xy^ for several values of 77 with S 23 = 0.4. 


R23{0,6) 


1 0 0 
0 cos 6 sin 0 
0 — sin 9 cos 9 


(3.51) 
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the PMNS matrix U in vacuum can be parametrized as 


u = /223(023,O) Rn{ei^, 5 ) i?12(012,O) 


C12C13 


S12C13 


si3e 




„i(5 


i<5 


'S12C23 — Cl2Sl3S23e C12C23 — Si2Sl3S23e C13S23 


„i(5 


iS 


'Sl2'S23 — Cl2Sl3C23e —C12S23 “ •Sl2'Sl3C23e C13C23 


(3.52) 


In the following, we rewrite the mixing matrix in matter U into the analogous form 

u = .^23(^23) 0)i?i3(6<i3, (5)iii2(0l2, 0) , (3.53) 

absorbing the extra mixing angles and CP phase into appropriate definitions of the ‘running’ 
parameters 9 12 , 0 i 3 , ^ 23 , and 6 . Frequent use is made of the relations 


Ri2{9,6)Q3 = Q3Ri2{9,5) , 
Ri3{0,m^ = ^ 3 ^ 13 ( 0 , 0 ) , 

^ 23 ( 0 , 0)g3 = Q3R23{0,-6) , 


(3.54) 


where Q 3 was defined in Eq. (3.4). 

• > 0 Case: 

Using Eq. (3.54), it is straightforward to show that 
U = UQ 3 VWX 

= R23{923, 0)Rl3i9l3, S)Ri2{9i2,0) Q 3 Rl2{<^, 0) R23{4>-, 0) .Rl2(X) “<5) 

^ ‘V' ^ ■v* ^ “V* ^ 

V V W X 

= R23{(^23, 0)Q3Ri3{6i3, O)i?i2(012, 0)i?i2(<^, 0)R23{4>7 0)Rl2{x, “<5) 

= R23{(^23, 0 )Q 3 Ri 3 { 6 i 3 , 0)Rl2i 9l2 + 0)i?23(<(', 0).Rl2(X) “ 

= ^12 

= -^ 23(^*235 O)g 3 i?i 3 ( 0 i 3 , O)i?i 2 ( 0 i 2 , 0 )i? 23 (</>) 0 )i 2 i 2 (X) —<^) > 


(3.55) 


where in the last and penultimate lines we have combined the two 12-rotations into 
one. We now commute R23{4>,0)Ri2{Xj^) through the other mixing matrices to the 
left as follows: 

— Step 1 : Commutation of R23{(p,0) through Ri 2 { 9 [ 2 , 0 ). 

In the range /? > 0, the angle 0(2 is approximately equal to tt/2 so we can 
approximate 

’’ 0 1 0’ 

.^12(^12)0) ~ .^ 12 (^ 2 ’*^) ^ 


-10 0 
0 0 1 


(3.56) 


Note that 


-R12 0) i223(</>, 0) — -Ri 3(</>, 0)i2i2 (^-,0^ 


(3.57) 
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for any (p. On the other hand, in the range /3 < — 1 the angle (p is negligibly 


small so we 

can approximate 





R2z{<Pi^) ~ .^ 23 ( 0 , 0 ) = 

'l 0 o’ 
0 10 

0 0 1 

= Ri3{0,0) . 

(3.58) 

Note that 

-^ 12(^125 0)-^23(0; 0) — .^ 13 ( 0 , O)i?l2(0i2) 0) 

(3.59) 

for any 0 '^ 2 - 

Therefore, for all (5 we have 





-Ri2(^'i 2! 0)-^23(</’) 0) PS i?i3((/>, O)i2i2(0i2,0) , 

(3.60) 

and 





U 

= .^23(^23) O)(53i?i3(0i3, d)Ri2{6'i2, 0)i?23(<(>) 0)iil2(X) —^) 

PS i?23(^23) O)(53i?i3(0i3, d)Ri2,{p, O)i?i2(0i2, 0)iii2(X) —^) 

= R23{923, 0)Q3Ri3{9i3 + <?), 0)i2i2(6'i2) 0)i?i2(X) “<^) 



= ^13 

= R23{023,0)Q5Ri5{e[s,0)Ri2{e[2,0)Rl2{X,-S) ■ 

(3.61) 


— Step 2 : Commutation of i?i2(X) ~^) through Ri2{0[2, 0). 

In the range /3 > 0, the angle 9^2 is approximately equal to 7r/2 as we have noted 
above and we have the approximation given in Eq. (3.56). Note that 

-^12 (-,o) iii2(x>-<5) = (^-,0^ (3.62) 

for any x- On the other hand, in the range /3 < 0, the angle x is negligibly small 
so we can approximate 


Ri2 {x,-S) ~ Rl2{0,-d) 


1 0 0 
0 1 0 
0 0 1 


.^ 12 ( 0 , 5 ) = i?23(0, —(5). (3.63) 


Note that 

Ri2{e[2,0)Ri2{0,-6) = i?12(O,5)i?l2(0'i2,O) (3.64) 

for any 9 'i2- Therefore, for all /3 we see that 

^12(6''i2>0)^12(x,-<^) ~ Ri2ix,S)Ri2i9u,0) , (3-65) 

and 

U « R23{923, 0)Q3RMs, 0)Ri2{X, 5)i?12(0'l2, 0) . (3.66) 
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— Step 3 : Commutation of Ri2{x-:^) through i?i 3 ( 0 '^ 3 ,0). 

When > 0 we have 0^3 ~ ^ in the range /3 > 1 so we can approximate 


-Ri3(<9'i 3,0) « i?i3(^-,0^ 


0 0 1 

0 1 0 

-10 0 


(3.67) 


Note that 

-^13 (-^, 0 ) i?i2(x,<^) = .R23 (X)-<^)-Ri3 (^ 2 ’^) (3.68) 

for any x- On the other hand, in the range /? < 0 the angle x was negligibly 
small so that we had Eq. (3.63). Note that 


-^13(^13) 0)-^12(0, 5) — -^ 23 ( 0 , —5)7?13(013, 0) , (3.69) 

for any Therefore, for all (3 we see that 

i?i3(0'i3>O)7?i2(x,5) « i223(x,-<5)72i3(0'i3>O) , (3.70) 

and using Eq. (3.54) we obtain 

U ^ R23{023,O)Q3R23{x,-S)RM3,O)Ri2{9[2,O) 

= R23{923, 0)i?23(X, O)7?13(0 'i3, ^)^12(0'i2,0)Q3 
= R23{923^, O)i?13(0'l3> '5)i?12(0'i2, 0)^3 

= 023 

= R23{023,O)Ri5{0[s,S)Rl2{0'i2,O)Q3 , (3.71) 


where in the last and penultimate lines we have combined the two 23-rotations 
into one. The matrix on the far right can be absorbed into the redefinitions 
of Majorana phases and can be dropped. 


Thus, we find that the effective mixing matrix U in the case Sin'll > 0 can be 
expressed as Eq. (3.53) with the effective mixing angles and effective CP-violating 
phase given approximately by 


12 " 

^0'l2 

= 012 

13 " 

^0'l3 

= 013 

23 " 

^ 023 

= 023 

6 ? 

s <5 . 



+ 7^ ) 

+ 4> ) 
+ x, 


(3.72) 


• < 0 Case: 


Using Eq. (3.54), we obtain 


U = UQsVWY 
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= -^23(^23, 0)-Ri3(^13, <5)-Ri2(^12, 0 ) Q 3 Rl2{‘f, 0 ) R23{4’^ 0 ) i?i3(—V’, S) 

' -V-' '-V--V--V-" 

U V W Y 

= R23(^23,O)Q3Ri3(0i3,O)Ri2(^l2>O)R23(4>jO)Ri3(—'!p,—S) . (3.73) 


We now commute R23((f>, 0)Ri3(—ip,—S) through the other mixing matrices to the 
left and re-express U as in Eq. (3.53), absorbing the extra mixing angles and CP 

I'N..' 

phase into 612, O13, 623, and 6. The first step is the same as the Sin'll > 0 case, the 
only difference being the /3-dependence of 0(3, which is also shown in Fig. 3(b). 


— Step 2 : Commutation of Ri3{—'ijj, —6) through Ri2{0[2, 0). 

In the range /3 > 0 the angle 0^2 is approximately equal to 7r/2 as we have noted 
previously, and we have the approximation given in Eq. (3.56). Note that 

Ri 2 (|, 0) Ri 3{-^, -d) = i?23(V’, -mi2 (|, 0) (3.74) 

for any V’- On the other hand, in the range /3 < 0 the angle is negligibly small 
so that 


i2i3(-V^,-5) « Ri3{0,-6) 


1 0 0 
0 10 
0 0 1 


R23{0,-6) = Ri3{0,6) . (3.75) 


Note that 

Ri2{0[2,0)Ri3{0,-S) = R23{0 ,-S)Ri2{0[2,0) . (3.76) 

Therefore, for all /3 we see that 

Ri2{0[2,0)Ri 3{-^P,-S) ^ R23ii^,-S)Ri2{0[2,0) , (3.77) 

and 

u ^ R23{e23,0)Q3Rl3{e[3, 0)R23ii^, -S)Rl2{0[2, 0) . (3.78) 


— Step 3 : Commutation of Ri3{ip, —5) through i?i3(033,O). 
When 6171^1 < 0 we have 0^3 ~ 0 in the range /3 > 1 so that 


-^13(^13)0) « i?i3(0,0) 


10 0 
0 1 0 
0 0 1 


(3.79) 


Note that 

-^13(0, O)i?23('0) —^) = 7?23(V'; —^).Ri3(0, 0) (3.80) 

for all 'll). On the other hand, in the range /? < 0 the angle if) was negligibly 
small so that we had the approximation Eq. (3.75). Note that 


^13(0'i 3, 0)7223(0,-5) = 7223(0,-<i)i2i3(0'i3>O) (3.81) 
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(a) Normal Hierarchy 


(b) Inverted Hierarchy 


Figure 8 . The /3-dependence of 6*23 for the (a) normal and (b) inverted hierarchies for several 
values of 77 with 5^3 = 0.4. 


for all O'l^. Therefore, for all /3 we see that 


^13(0'l3>O)i?23(V',-<5) « i?23(V’,-<5)^13(0'l3,O) , (3-82) 

and using Eq. (3.54) we obtain 

U « /?23(023,O)Q3i223(V’,-<5)i2l3(0'l3:O)i?12(0'l2,O) 

= ^23(023, 0)i?23(V', O)iil3(0'l3, mi2{0'i2. 0)^3 
= ^23(023j^, O)i?13(0'l3> 5)i?l2(0'i2, 0)^3 

= 023 

= R23i023,0)Rni0[3, S)Ri2i0[2,0)Q3 , (3.83) 


where in the last and penultimate lines we have combined the two 23-rotations 
into one. The matrix Q 3 on the far right can be absorbed into redefinitions of 
the Majorana phases and can be dropped. 

Thus, we find that the effective mixing matrix U in the case 6 m^i < 0 can be 
expressed as Eq. (3.53) with the effective mixing angles and effective CP-violating 
phase given approximately by 


012 ~ 012 ~ 012 + f ) 

013 « 0'i3 = 9i5 + 4> , 

023 « 023 = 023 + V’ 5 

5 « 5 . 


(3.84) 


3.4 Summary of Neutrino Case 

To summarize what we have learned, inclusion of the term in the effective Hamilto¬ 

nian shifts 023 to 023 = 023 + X tti® 5m\i > 0 case, and to 023 = 023 + '*/’ for the Smh < 0 
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case. For both cases, ^23 can be calculated directly without calculating x or V' first via the 
expression 


tan 20 n 


23 


[^"^31^13 - <^"^2i(ci2 - gL’Sp] sin 26^23 

cos2023 - ‘^arj 


(3.85) 


TT 


Note that as /3 is increased, 023 runs toward — if Sm^iT] > 0, while it runs toward 0 if 
< 0. The /3-dependence of 0^3 is shown in Fig. 8 . The CP-violating phase 6 is 
unaltered and maintains its vacuum value. 

The running of the effective mass-squared differences are also modified in the range 
/3 > 0. For the > 0 case, Ai and A 2 show extra running, while for the < 0 case, 
it is Ai and A 3 that show extra running, cf. Figs. 6 and 7. 


3.5 Discussion at the Probability Level 

So far we have focused our attention on how the flavor-diagonal NSI parameter rj modifies 
the running of the effective mass-squared differences, mixing angles, and CP-violating phase 
as functions of a = a(l -|- C) for the neutrinos. We derived simple analytical expressions 
for these effective parameters using the Jacobi method. We have discussed the running of 
these effective neutrino oscillation parameters for both the normal {5m\i > 0 ) and inverted 
< 0) neutrino mass hierarchies. The modifications induced by rj and C the run¬ 
ning of effective oscillation parameters in the case of anti-neutrino are discussed in detail 
in appendix A. At this point, we look at how these lepton-flavor-conserving NSI parame¬ 
ters alter the neutrino oscillation probabilities for various appearance and disappearance 
channels. 

In the three-flavor scenario, the neutrino oscillation probabilities in vacuum^ for the 
disappearance channel (initial and final flavors are same) take the form 

P{Ua -)■ Ua) = I- 4 |t/„2p (l “ \Ua2\‘^) sin^ ^ “ 4 |C/„3p (l - \Ua3\‘^) sin^ ^ 
-F2|t/„2p|t7a3p ^4sin^ ^sin^ ^-ksinA2isinA3i^ , (3 

and for the appearance channel (initial and final flavors are different) we have 

P{i'a = 4 1142Pi 17/32P sin^ ^ -h 4 |t/„3p|I7^3p sin^ ^ 

sin A 31 
n A 2 i^ 



+2 J? {U*^Ui33Ua2U^2) ^4sin^ sin^ -|- sin A 21 
-43= {Ul^UP 2 .Ua 2 Up 2 ) ^sin^ ^ sin A 31 - sin^ ^ si 



In the above equations, we define 


Aii = —PL = 2.534 
^ 2E 




tA(GeY\f L \ 


6 m% = mf — mi 


eV^ J \ E J \ km J 




(3.88) 


®We follow the conventions, notations, and the expressions for various neutrino oscillation probabilities 
in vacuum as given in the appendices A.l and A.2 of Ref. [70]. 
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The transition probabilities in Eq. (3.86) and Eq. (3.87) contain several elements of the 
PMNS matrix U which are expressed in terms of three mixing angles 0i2, ^ 23 ) and a CP- 

violating phase 6 as shown in Eq. (3.52). The oscillation probabilities for the anti-neutrinos 
are obtained by replacing Uai with its complex conjugate. 

The neutrino oscillation probabilities in the presence of matter are obtained by replac¬ 
ing the vacuum expressions of the elements of the mixing matrix U and the mass-square 
differences Ajj with their effective ‘running’ values in matter [21-23] : 

Uai Uai (^12 ^12)^13 ^13)^23 ^ 23 ) ) ^ij ^ij ~ —’ (3-89) 

and for the anti-neutrinos 

Uai Uai {012 ^125^13 ^13)^23 ^ 23 ) ) ^ij ^ij ~ —' (3-90) 

To demonstrate the accuracy (or lack thereof in special cases) of our approximate analytical 
results, we compare the oscillation probabilities calculated with our approximate effective 
running mixing angles and mass-squared differences with those calculated numerically for 
the same baseline and line-averaged constant matter density along it. Eor the mixing angles 
and mass-squared differences in vacuum, we use the benchmark values from Ref. [194] 
as listed in Table 1. In some plots, we take different values of sin^ 023 and 6 which we 
mention explicitly in the figure legends and captions. In this paper, all the plots (except 
in appendix B) are generated considering the line-averaged constant matter density for 
a given baseline which has been estimated using the Preliminary Reference Earth Model 
(PREM) [195]. In appendix B, we compare the exact numerical probabilities with line- 
averaged constant Earth density and varying Earth density profile for 8770 km and 10000 
km baselines. 

In Eig. 9, we present our approximate —)■ t'e oscillation probabilities (blue curves) as 
a function of the neutrino energy against the exact numerical results (red curves) consider¬ 
ing® r] = 0.1, C = 0 (left panels) and r] = —0.1, C = 0 (right panels). The upper panels are 
drawn for the baseline of L = 2300 km, which corresponds to the distance between CERN 
and Pyhasalmi [196-198] with the line-averaged constant matter density of p = 3.54 g/cm^. 
In the lower panels, we give the probabilities for the baseline of L = 8770 km, which is the 
distance from CERN to Kamioka [199] assuming p = 4.33g/cm^. Here, in all the panels, 
we assume sin^ 023 = 0.41 (023 = 40°), 5 = 0°, and normal mass hierarchy { 6 m\i >0). To 
see the differences in the oscillation probability caused by the NSI parameters, we also give 
the exact numerical SM three-flavor oscillation probabilities in matter in the absence of 
NSTs which are depicted by the solid black curves with figure legend ‘SM, Exact.’ It has 
been already shown in Ref. [70] that our approximate expressions for the p = 0 case match 
extremely well with the exact numerical results for all these baselines and energies. We 
also compare our results with the approximate expressions of Asano and Minakata^ [69] 

®We take (^ = 0 in our plots since we expect it to be hidden in the uncertainties in the matter density 
and neutrino energy. 

^For comparison, we take Eq. (36) of Ref. [69] where the authors adopted the perturbation method 
[80, 200] to obtain the analytical expressions for oscillation probability in presence of NSI’s for large 6 ^ 13 . 
The same analytical expressions are given in a more detailed fashion in Eqs. ( 8 ) to (13) in Ref. [201]. 
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L=8770 km, (£^^=0.10, €„=-0.l0)/(n=0.l0, (= 0) L=8770 km, (e^^=-0.10, e„=0.10)/(;;=-0.10, ^=0) 




Figure 9. ^ Ve transition probability as a function of neutrino energy E in GeV for 2300 km 

(8770 km) baseline in upper (lower) panels. We compare the analytical expressions of this work and 
Asano-Minakata [69] against the exact numerical result assuming ry = 0.1, C = 0 (left panels) and 
ry = —0.1,^ = 0 (right panels). The solid black curves portray the standard three-flavor oscillation 
probabilities in matter without NSI’s. In all the panels, we consider 6*23 = 40°, (5 = 0°, and normal 
mass hierarchy. 


(dashed green curves). The correspondence between our rj and ^ and the NSI parameters 
See) S/i/i) and Err used in Refs. [69, 201] can be obtained via Eqs. (2.4), (2.7), and (2.9) 
which suggest the changes : a —)■ a = a(l -|- C}, See —> 0, —)• ry, and Err —^ —f]- We can 

see from Fig. 9 that for the 2300 km baseline, the Asano-Minakata expressions give better 
matches compared to our results, while for the 8770 km baseline, our expressions agree 
better with the exact numerical results. 

The accuracy of our analytical approximations as compared to the exact numerical 
results for different vacuum values of 623 is demonstrated in Fig. 10. Here we consider the 
minimum (35°) and maximum (55°) values of 823 which are allowed in the 3a range [194]. 
We also present the results for the maximal mixing choice. All the plots in Fig. 10 have 
been generated assuming (5 = 0° and normal mass hierarchy { 6 m\i > 0). We consider the 
same choices of i] and C as in Fig. 9 and results are given for 2300 km (upper panels) and 
8770 km (lower panels) baselines. As is evident, our approximation provides satisfactory 
match with exact numerical results for different values of 023 - 


- 29 - 








L=2300 km, (e.^=0.10, e^^=-0.10)/(77=0.10, 0) 


L=2300 km, (e^^=-0.10, err=0.10)/(?7=-0.10, (= 0) 



E (GeV) 

L=8770 km, {e^^=0.10, e^^=-0.10)/{7/=0.10, C= 0) 


35°, Exact 

— ^ 23 = 35°. This work 
__ _ ^ 23 = 45°, Exact 
_ 023 = 45°. This work 
; _ 02.3= 55°, Exact 

_ 023= 55°, This work 


E (GeV) 

L=8770 km, (e^p=-0.10, e,,=0.10)/(7/=-0.10, f= 0) 



Figure 10. Comparison of our analytical expressions (solid curves) to the exact numerical results 
(dashed curves) for various values of 023 assuming 5 = 0° and 5m\i > 0. Upper (lower) panels are 
for 2300 km (8770 km) baseline. 


L=2300 km, (e^^=0.10, e„=-0.10)/(;7=0.10, f= 0) L=2300 km, (e^^=-0.10, e„=0.10)/(;;=-0.10, f= 0) 



E (GeV) E (GeV) 

Figure 11. Comparison of our analytical expressions (solid curves) to the exact numerical results 
(dashed curves) for four different values of the CP-violating phase 5 at 2300 km assuming 023 = 40° 
and > 0. 

Fig. 11 compares our approximate probability expressions (solid curves) against the 
exact numerical results (dashed curves) assuming four different values of the CP-violating 
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L=8770 km, (e^^=0.10, 6^^=-0.10)/(7/=0.10, (= 0) L=8770 km, (e^^=-0.10, €^^=0.10)/(?7=-0.10, 0) 




(a) 77 = 0.1,C = 0 (b) 77 =-0.1,C = 0 

Figure 12. ^ survival probability as a function of neutrino energy E in GeV for two 

different values of 023 at 8770 km baseline. Comparison between the analytical and numerical 
results assuming 77 = 0.1, ^ = 0 (left panel) and 77 = — 0.1,C = 0 (right panel). The standard 
three-flavor oscillation probabilities in matter without NSI’s are also shown. In both the panels, we 
assume (5 = 0° and Smh > 0. 

phase 6 (0, 90, 180, and 270 degrees) at 2300 km. Here, we consider 623 = 40° and > 0. 
These plots clearly show that our approximate expressions work quite well even for non¬ 
zero 6 and can predict almost accurate LjE patterns of the oscillation probability for finite 
5 and 77 . It also suggests that one can explain qualitatively the possible correlations and 
degeneracies between <5 and r/ using these analytical expressions which cannot be tackled 
with numerical studies. 

In Fig. 12, we plot the —>■ survival probability in the presence of NSI for two dif¬ 
ferent vacuum values of O23 (40° and 50°) at 8770 km. We show the matching between the 
analytical and numerical results assuming rj = 0.1, C = 0 (left panel) and rj = —0.1,C = 0 
(right panel). We also give the exact numerical standard three-flavor oscillation proba¬ 
bilities in matter without NSI’s so that one can compare them with the finite 77 case. In 
both the panels, we assume 5 = 0° and > 0. Fig. 12 shows that our approximate 

expressions match quite nicely with the numerical results. Note that at higher energies, the 
impact of NSFs are quite large in the —)• survival channel and there is a substantial 
difference in the standard and NSI probabilities for both the choices of 623 which can be 
probed in future long-baseline [202, 203] and atmospheric [204-206] neutrino oscillation 
experiments. 

Another important point to be noted that in the absence of NSI, the standard proba¬ 
bility curves for both the values of 623 almost overlap with each other at higher energies, 
whereas with NSI, there is a large separation between them. It immediately suggests that 
the corrections in the probability expressions due to the NSI terms depend significantly 
on whether the vacuum value of 623 lies below or above 45° [31]. Fig. 12 also indicates 
that there are degeneracies between the octant of 623 and the sign of NSI parameter rj for 
a given choice of hierarchy. For an example, P^^(023 = 50°, 77 = 0.1) in the left panel 
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(a) Normal and Inverted Hierarchy 


(b) Normal Hierarchy 


Figure 13. /3-dependence of (a) sin^ ( 205 . 3 ) the normal {Sm^i > 0) and inverted {Srn^i < 0) 
hierarchies, and (b) that of 5^3 for various values of ij for the normal hierarchy case. The dashed 
peak in (b) indicates the behavior of sin^( 2053 ) when Srri^i > 0 . 


is almost same with P^^{023 = 40°, r/ = —0.1) in the right panel for the energies above 
12 GeV or so. Again, P^^{023 = 40°, r/ = 0.1) in the left panel matches quite well with 
= 50°, T/ = —0.1) in the right panel. These kinds of degeneracies can be well ex¬ 
plained qualitatively with the help of our analytical expressions. We discuss this issue in 
detail in the next section which is one of the highlights of this work. 

4 Possible Applications of Analytical Expressions 

In this section, we discuss the utility of our analytical probability expressions to determine 
the conditions for which the impact of the NSI parameter rj becomes significant. We also 
give simple and compact analytical expressions to show the possible correlations and de¬ 
generacies between 623 and rj under such situations. We begin our discussion with electron 
neutrinos. 

4.1 r'e —t'a Oscillation Channels 

Let us hrst consider the I'e —)■ i^a (a = e, //, r) oscillation channels in matter in the presence 
of the NSI parameter rj. Since we expect the effect of r] to become important in the range 
/3 > 0, we set 5^2 ~ 1, c '12 ~ 0 (which is valid in the range j3 S> —2, see Fig. 2(a)), which 
leads to the following simple expressions: 




P{T^e 

—)■ Vf.) - 

si — sin^(2053) 

• 2 ^32 

sm - , 

2 ’ 

(4.1) 



s S23 sin^(20i3) 

. 2 ^32 
sm 2 ’ 

(4.2) 


w) - 

s C23 sin^(2053) 

. 2 ^32 

sm - . 

2 

(4.3) 
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L=^OQOQ km, 023=40°, /7=0.1, Normal Hierarchy L=^0000 km, 023=40°, /7=0.1, Normal Hierarchy 



(a) (b) 


Figure 14. The oscillation probabilities (a) P{i'e —t and (b) P{i'e —t for the normal 

hierarchy case with = 0.4. The red lines are the standard oscillation probabilities with rj = 0 
and 6 = 0,90,180, and 270 degrees. The blue lines are our approximate analytical results with 
77 = 0.1 and the black dashed lines are the numerically calculated probabilities with 77 = 0.1. In 
both the cases, we plot the probabilities for four values of i5 = 0,90,180, and 270 degrees. Note 
that P{iye —t t'/j) is enhanced by about 77 , while P{ue —>■ Vr) is suppressed by an equal amount. 

Recall that when > 0, the effective mixing angle increases monotonically toward 
7r/2, going through vr/d around /3 ~ 0, while in the Smh < 0 case, it decreases mono¬ 
tonically toward 0, cf. Fig. 3(b). This will cause sin^( 20 g 3 ) to peak prominently around 
/3 ~ 0 for the normal hierarchy case, but not for the inverted hierarchy case as shown 
in Fig. 13(a). As discussed in Ref. [70], demanding that sin^(A32/2) also peaks at the 
same energy leads to the requirements of L ~ 10000 km and E ^ 7 GeV. Thus, measuring 
P{i'e t'e) survival probability at this baseline and energy will allow us to discriminate 
between the normal and inverted mass hierarchies irrespective of the value of 023 [207]. 
Also, around /3 ~ 0, A 32 is not affected by rj provided 77 < 0.1 (see middle and bottom 
panels of Fig. 6 ), allowing this channel to determine the mass hierarchy free from any NSI 
effect. 

To see the effect of 77 , we need to observe the running of 623 - This could be visible in 
the I'e —^ t'/i and Ve appearance channels for the normal hierarchy case as a change 

in the heights of the oscillation peaks around E ^ 7 GeV provided 023 deviates sufficiently 
from the vacuum value of 023 at that energy. The running of S 23 for various values of 77 is 
depicted in Fig. 13(b). At d ~ Eq. (3.85) tells us that 

s'is ~ ■523 + ^) ( 4 - 4 ) 

showing the possible shift in 023 due to 77 in a simple and compact fashion which clearly 
establishes the merit of our analytical approximation. Eq. (4.4) also shows compactly pos¬ 
sible correlations and degeneracies between 623 and 77 . Such correlations and degeneracies 
could also be found numerically, but the reason for those features will not be so transpar¬ 
ent. Note that Eq. (4.4) suggests that a positive value of rj would enhance P{ire —>• ir^) 
while suppressing P{i'e and a negative rj would do the opposite. Eig. 14 confirms 
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this feature where we plot the Vf. —)• {Vf. —>• v^) transition probability in the left (right) 
panel. In both panels, the standard oscillation probabilities in matter without NSI’s for 
four different values of CP phase® S are given by the solid red lines. Blue lines (black dashed 
lines) depict the approximate analytical (exact numerical) probabilities with rj = 0.1 and 
four different values of 5. Fig. 14 infers that the effect of rj is visible in these channels 
provided the vacuum value of 623 is sufficiently well known and rj is also large enough. 

4.2 —)> Va Oscillation Channels 

Next, let us consider —)■ t'a (a = /r, r) oscillation channels. In addition to assuming 
(3 S> — 2 , which allows us to set s '^2 ~ 1) c '^2 ~ 0) further restrict our consideration to 
the range /3 > 1 , which allows us to set s '^3 1 , c '^3 0 or 5)^3 0 , 0)^3 « 1 depending 

on whether > 0 or < 0 (see Fig. 3(b)). With these conditions, we obtain the 

following simple expressions: 


- S13S2I) sin^ ^ - 4 c;| 4 |(l - cf^sg) sin^ ^ 


+2'Si3Ci 3'S23 ( 4 sin" 


Ao 






1 — sim( 2023 ) sin- 
1 — sin^( 2023 ) sin' 


2 A31 


2 A 31 ~ ~ 1 

sin 2 -h sin A 21 sin A 31 1 

(4.5) 

( 6 ml. > 0 ) 

< 0 ) 

(4.6) 


P{v^ -t t^t) ~ sin( 26 » 23 ) 


/4 . 2 A2I , /4 ■ 2 ^31 

•Sis sin — + Ci 3 sm — 


„/2 /2 






sin^ ( 26 * 23 ) sin" 
sin^ (26*23) sin" 


2 sin^ 

A21 0 A31 1 ''y T \ 

2 sm^ ^ + 2 1 

(4.7) 

A21 

{ 6 ml^ > 0) 

(4.8) 

^2 

A31 

2 

{ 6 ml^ < 0) 



In the absence of rj, 023 does not run and sin^(2023) will maintain its vacuum value close to 

TT 

one. In the presence of a non-zero rj, however, 023 will run towards either — or 0 depending 
on the sign of 6 m\iri, as was shown in Fig. 8 , and sin^(2023) will run toward zero in both 
cases. This is depicted in Fig. 15. 


®Fig. 14 shows that the impact of the CP phase <5 is quite weak around /3 ~ 0. In this region, 9'^2 
approaches to 7r/2 so that s '12 « 1 and C 12 ~ 0. Therefore, the Jarlskog Invariant [208] in matter J = 
Si 2 Ci 2 Si 3 ci 3 S 23 C 23 sin5 almost approaches to zero diminishing the effect of 5. This argument works even in 
the presence of the NSI parameter rj. Our simple and compact approximate probability expressions given 
by Eq. (4.2) and Eq. (4.3) also validate this point as there are no 5-dependent terms in these expressions. 
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p p 

(a) Normal Hierarchy (b) Inverted Hierarchy 

Figure 15. The running of sin^(2023) foi' various values of 77 for the (a) normal and (b) inverted 
mass hierarchies with s^g = 0.4. The blue/green lines indicate positive/negative values of rj. 



P p 


(a) 


(b) 


Figure 16. The running of (a) A 21/2 and (b) sin^(A 2 i/ 2 ) for various values of rj for the normal hier¬ 
archy case with L = 10,000 km and p = 4.53 g/cm^. The blue/green lines indicate positive/negative 
values of 77 , while the red lines are drawn for 77 = 0 case. 


Let us see how the A factors in Eq. (4.6) and Eq. (4.8) behave in the range /3 ^ 1- For 
the normal hierarchy case, we have 


(5A21 


Ax+ — ^"x- ~ 2 a|r 7 | (77 7 ^ 0 ) 

\"_ - A'_ PS (77 = 0) 


while for the inverted hierarchy case, they take the form 


(5A31 


Ay_ - \y+ ~ - 2 a 1 77 1 (^ / 0 ) 

- A'_ ^ -\5ml^\cl^ {r] = 0) 


(4.9) 


(4.10) 


Since the sign of 6 Xij does not affect the value of sin^(Aij/ 2 ), both mass hierarchies lead to 
the same asymptotic oscillation probabilities. We will therefore only consider the normal 


- 35 - 






















023 ©23 



Figure 17. 023 sin^(2023) as a function of positive and negative values of r] for a fixed S 23 = 0.4 

(sis = 0.6) in the left (right) panel. The upper (lower) signs are applicable for the normal (inverted) 
hierarchy case. 


hierarchy case in the following. Recalling (see Eq. (2.2)) that 


a = 2V2GFNeE = 7.6324 x 10“^(eV2) 
(where we have set C = 0) and 


Aij = 2.534 


= 1.934 


dA. 




eV" 


<5 A, 


V 


we find 
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<^A 21 
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10^ km 
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g/cm" 


(¥)(a 


L 


10^ km 


30 


L 




\ E J V^lO'^km 


(r/ / 0) 

iv = 0 ) 


(4.11) 


(4.12) 


(4.13) 


Therefore, for fixed baseline L and matter density p, as the neutrino energy E is increased, 
A 21/2 damps to zero when rj = 0, but asymptotes to a constant value proportional to \r]\ 
when r] ^ 0. This is demonstrated in Fig. 16 for the baseline L = 10000 km with average 
matter density p = 4.53g/cm^. Consequently, when r/ = 0, the factor sin^(2023) stays 
constant while sin^(A2i/2) damps to zero as we increase E, while in the p ^ 0 case, the 
factor sin^(2023) damps to zero while sin^(A2i/2) asymptotes to a constant value as E is 
increased. In either case, the —)• u-j- oscillation probability is suppressed at high energy. 
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The difference between the rj = 0 and rj ^ 0 cases could manifest itself at the —?• Ur 

oscillation peak which happens at 


TT 

2 





20GeV , 15 


0.6 , 


(4.14) 


for L = 10000 km, p = 4.53 g/cm^ case. (5 ~ 0.6 is on the borderline of the applicability of 
our ,0 > 1 approximation. Nevertheless, let us examine what our approximation suggests. 
Expanding Eq. (3.85) for small ap, we find 


^23 ~ ^23 + 


asin(2023) 




p + 


13 


(4.15) 


which at j5 ~ 0.6 yields 


0^3 023 ± 3r? . (4.16) 

In the above equation, upper (lower) sign is applicable for the normal (inverted) hierarchy 
case. It suggests that there is a degeneracy between the choices of sign of and the 
sign of p which give rise to same amount of corrections in 023- To observe the shift in the 
oscillation probability, we keep terms up to order p'^ since the linear term is suppressed 
due to the fact that 023 is close to 7 r/ 4 . Therefore, we have 


sin(2023) = sin(2023) + 2 cos(2023)(5023 - 2 sin(2023)((5023)^ H- 

« sin(2023) ± 6cos(2023)r/ - 18 sin(2023)7?^ . (4.17) 

Eor the benchmark value of 5^3 = 0.4 (0.6) in the lower (higher) octant, above equation 
can be written as 


sin(2023) 


sin(2023) ±p- I8p‘^ ( 5^3 = 0.4) 

sin(2023) T 0 - 180 ^ (sig = 0.6) ’ 


(4.18) 


where upper (lower) signs are for the normal (inverted) hierarchy. The above equations 
clearly reveal that for a given choice of hierarchy, there are degenerate solutions® of the 
octant 0/623 and the sign of NS I parameter p (S 23 = 0 . 4 , 7 ] = ±0.1 and S 23 = 0.6, p = ±0.1), 
giving rise to same value of sin(2023). Eig. 17 shows the variation in 0^3 and sin^(2023) 
as a function of positive and negative values of p for a fixed ~ ('^23 ~ 0 . 6 ) in the 

left (right) panel. The upper (lower) signs correspond to the normal (inverted) hierarchy 
scenario. We can see from the left panel of Fig. 17 that assuming normal hierarchy and 
S 23 = 0.4, the value of sin^(2023) gets reduced by substantial amount for p = —0.1 case as 
compared to p = 0.1 as given by Eq. (4.18). It means that negative values of p, which shifts 
023 further away from 7 r/ 4 , would lead to a larger suppression of the —)■ Ur oscillation 
probability, and enhancement of the —)• survival probability. The left and right panels 
of Fig. 18 exactly show this behaviour where we plot the approximate analytical and exact 
numerical —?• Vr (left panel) and —?■ (right panel) probabilities for p = ± 0.1 

assuming S 23 = 0.4, > 0 and L = 10000km, p = 4.53g/cm^. The situation gets 

reversed completely for cases in which S 23 > 0-5 which is quite evident from the right panel 
of Fig. 17 where we consider = 0.6. All these observations in Fig. 17 and Fig. 18 suggest 
that our approximate calculations are valid. 

®We have already seen this degeneracy in the —>■ t'^i) oscillation channel in Fig. 12. 
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Figure 18. (a) —?> v-y and (b) —>■ oscillation probabilities for s^g = 0.4, L = 10,000 km, 
and p = 4.53 g/cm^. The red lines are the standard oscillation probabilities with ry = 0. The solid 
blue (green) lines are our approximate analytical results with 77 = 0.1 (77 = —0.1). The black dashed 
lines are the numerically calculated probabilities. In both the panels, we take (5 = 0° and (5m|g > 0. 


5 Summary and Conclusions 

Analytical studies of the neutrino oscillation probabilities are inevitable to understand how 
neutrino interactions with matter modify the mixing angles and mass-squared differences in 
a complicated manner in a three-flavor framework. In previous papers [68, 70], we showed 
that the neutrino oscillation probabilities in matter can be well understood if we allow 
the mixing angles and mass-squared differences in the standard parametrization to ‘run’ 
with the matter effect parameter a = ‘l\f 2 GpN(.E^ where Nf. is the electron density in 
matter and E is the neutrino energy. We managed to derive simple and compact analytical 
approximations to these running parameters using the Jacobi method. We found that 
for large 613 , the entire matter effect could be absorbed into the running of the effective 
mass-squared differences and the effective mixing angles 612 and 613 , while neglecting the 
running of the mixing angle 623 and the CP-violating phase <5. 

In this paper, we extended our analysis to study how the running of the neutrino 
oscillation parameters in matter would be altered in the presence of NSI’s of neutrinos 
with the matter fermions. Such NSI’s are predicted in most of the new physics models 
that attempt to explain the non-zero neutrino masses, as well as in a wealth of various 
other BSM models. There, the NSI’s are simply the effective four-fermion interactions at 
the energy scales relevant for neutrino oscillation experiments that remain when the heavy 
mediator fields of the full theory are integrated out. These NSI’s give rise to new neutral- 
current type interactions of neutrinos (both flavor-conserving and flavor-violating) during 
their propagation through matter on top of the SM interactions, causing the change in the 
effective mass matrix for the neutrinos which ultimately affect the running of the oscillation 
parameters and hence change the oscillation probabilities between different neutrino flavors. 
These sub-leading new physics effects in the probability due to NSTs can be probed in 
upcoming long-baseline and atmospheric neutrino oscillation experiments. 
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In this work, we restricted our attention to the matter effect of flavor-conserving, non- 
universal NSI’s of the neutrino, relegating the discussion of the flavor-violating NSI case 
to a separate paper [73]. The relevant linear combinations of the flavor-diagonal NSI’s 
were rj = — £tt)I‘^ and C = ^ee — where a non-zero C led to a rescaling 

of the matter-effect parameter a —>• a = a(l-|-C)) while a non-zero rj led to non-trivial 
modifications on how the running oscillation parameters depend on a. 

Utilizing the Jacobi method, as in Refs. [ 68 , 70], we obtained approximate analytical 
expressions for the effective neutrino oscillation parameters to study how they ‘run’ with 
the rescaled matter-effect parameter a, and to explore the role of non-zero t] in neutrino 
oscillation. We found that in addition to the two rotations, which were required for the SM 
matter interaction and were absorbed into effective values of 0 i 2 and 613 , a third rotation 
was needed to capture the effects of rj, which could be absorbed into the effective value 
of 023- Thus, within the neutrino mixing matrix, the effect of rj appears as a shift in the 
effective mixing angle 823 , while the SM matter effects show up as shifts in 612 and ^ 13 . 
The CP-violating phase 5 remains unaffected and maintains its vacuum value. The running 
of all the effective neutrino oscillation parameters were presented for both the normal and 
inverted neutrino mass hierarchies. The changes caused by r] in the running of the effective 
oscillation parameters for the anti-neutrino case are discussed in detail in appendix A. 

We have also studied the impact of the lepton-flavor-conserving NSI parameters on the 
neutrino oscillation probabilities for various appearance and the disappearance channels. 
To demonstrate the accuracy (or lack thereof in special cases) of our approximate analyti¬ 
cal expressions, we compared the oscillation probabilities estimated with our approximate 
effective ‘running’ mixing angles and mass-squared differences with those calculated nu¬ 
merically for the same choices of benchmark oscillation parameters, energy, baseline, and 
line-averaged constant matter density along it. We found that our approximation provided 
satisfactory matches with exact numerical results in light of large ^13 for different values 
of O23, CP-violating phase 6, and for positive and negative values of the NSI parameter r]. 
A comparison of our results with the approximate expressions of Asano and Minakata [69] 
for the —)• zzg appearance channel has also been presented. 

Finally, we examined the merit of our analytical probability expressions to identify the 
situations at which the impact of the NSI’s become compelling. It was found that at higher 
baselines and energies, the impact of rj can be quite significant in the —)> survival 
channel if \r]\ is of the order of its current experimental upper bound. A considerable 
difference between the SM and NSI probabilities can be seen irrespective of the vacuum 
value of 023 , and the sign of the NSI parameter rj. We note that this feature may be 
explorable with the upcoming 50 kiloton magnetized iron calorimeter detector at the India- 
based Neutrino Observatory (INO), which aims to detect atmospheric neutrinos and anti¬ 
neutrinos separately over a wide range of energies and path lengths [209]. Using our 
analytical approach, we showed in a very simple and compact fashion that the corrections 
in 023 due to the rj depend significantly on whether the vacuum value of 023 lies below or 
above 45'’, suggesting a possible degeneracy between the octant of 023 and the sign of r] for 
a given choice of mass hierarchy. 
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A Effective Mixing Angles and Effective Mass-Squared Differences 
— Anti-Neutrino Case 


In this appendix, we study the matter effect due to the anti-neutrino NSI’s. We again 
utilize the Jacobi method to estimate how the NSI parameter r] alters the ‘running’ of the 
effective mixing angles, effective mass-squared differences, and the effective CP-violating 
phase 6 in matter for the anti-neutrinos. Like the neutrino case, we also present here a com¬ 
parison between our approximate analytical probability expressions and exact numerical 
calculations towards the end of this appendix. 


A.l Differences from the Neutrino Case 


For the anti-neutrinos, the effective Hamiltonian is given by 



1 

o 

o 

rH 

1_ 


1 

0 

0 

0 

_1 


’ 100 ’ 


1 

0 

0 

0 

1_ 

Hr, = U* 

0 A 2 0 

0 0 A 3 

u^ = u* 

0 0 

0 0 

-a 

0 0 0 

0 0 0 

—arj 

0 1 0 

0 0-1 


V V 

= Mg ^ =M^ 


= Ho 


(A.l) 

The differences from the neutrino case are the reversal of signs of the CP-violating phase 
6 (and thus the complex conjugation of the PMNS matrix U), and the matter interaction 
parameter a = a(l -|- C)- We denote the matter effect corrected diagonalization matrix as 
U (note the mirror image tilde on top) to distinguish it from that for the neutrinos. 


A.2 Diagonalization of the Effective Hamiltonian 
A.2.1 Change to the Mass Eigenbasis in Vacuum 

Using the matrix Qo from Eq. (3.4), we begin by partially diagonalize the effective Hamil¬ 
tonian Hn as 
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0 0 

-aQoU^ 

0 0 0 

U*Ql-ar,QoU^ 
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(MCO 
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U*Ql , (A. 2 ) 


M. 


= M^(6'i2,013,6'23) 


= Hn 


Mr, 


— M {9i2,Oio, 6 * 23 , ^) 
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M1(012,013,023) 



Ueiu:; 


^12^13 C12S12C13 C12C13S13 
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Cl2Cl3'Sl3 ■S12C13S13 5^3 
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— (6*12,6*13,6*23) , (A.3) 
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sin(26*i2) sin(26*23)si3 cos 5 + (^12 - C?2 Si 3 ) COs( 26*23) 
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12 '= 


— c 


12 


e *'^)si 3 sin(26*23) - (1 + Si3)'Si2Ci2 003(26*23) 


-S12C13 sm(26*23)e+*^ + C12S13C13 003(26*23) 
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— M^* (012, 013, 023, ■ 


(A.4) 


A.2.2 77 = 0 Case, First and Second Rotations 

As in the neutrino ease, we will first approximately diagonalize h'q and then add on the 
arjM'^ term later. The Jaeobi method applied to h'q is as follows; 

1. First Rotation 

Define the matrix V as 


V = 


C/^ 0 

0 0 1 


where 

— COS j 

Then, 


= sin If , tan 2 ip = — 


acf^ sin 26*12 


Sm 2 i + acfg eos 26*12 


(A.5) 


(- 6*12 < (/? < 0 ) . 

(A. 6 ) 


h” = v^h',v 
^ a'_ 

0 




-ac'^20l3'Sl3 

-aSi20l3'Sl3 


-ac'^20l3'Sl3 


■Q'Si2013'S13 


+ (5mq 


31 


where 


and 


c'i 2 = eos 6*12 , 


S 12 = sin 6*12 , 6*12 = 6*12 + (/? 


A+ = 


(Smh - acfg) ± C( 5 m^i + ac^s)^ - 400^35^2 <^"*-2 


21 


(A.7) 

(A. 8 ) 

(A.9) 
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Figure 19. (a) The dependence of 9^2 on ^ = — logj(a/|(5m|^|). (b) The /3-dependence of Aj 


The angle 9 'i 2 = ^12 + ^ can be calculated directly without calculating ip via 


tan 20^2 = 


sin 2612 
cos 2612 + ac 


13 


(0 < 9^2 < 9i2) 


(A.IO) 


The dependences of 9^2 and on /3 are plotted in Fig. 19. 

Note that in contrast to the neutrino case, 9^2 decreases monotonically from 9 i 2 to 
zero as /3 is increased. The /3-dependences of s '^2 = sin0^]^2 and 0^2 = cos0^]^2 are 
shown in Fig. 20(a). As /3 is increased beyond /3 = —2, C 12 grows quickly to one 
while S 12 damps quickly to zero. The product asi 2 stops increasing around /3 = —2 
and plateau’s to the asymptotic value of 617121812012 /= |(5m|^|0(e^) as shown in 
Fig. 20(b). That is: 


as'i2 = \6mli\0{e 2^) < \6ml^\0{e^) , 

ac'i2 = |<3mii|0(e"^) . (A.ll) 

Note also that the scales of A^^ are given simply by 

A'_ = -0(a) = 0(e"^) , 

A'+ = 0(,5mii) = |5mii|0(e2) , (A.12) 


since no level crossing occurs in this case. 


2. Second Rotation 


Since adi 2 continues to increase with /3 while as '^2 does not, we perform a (1,3) 
—n 

rotation on Rg next. 

Define the matrix W as 


W = 


Ocj) 0 Sfp 
0 1 0 


(A.13) 
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Figure 20. (a) The /3-dependence of s '12 = sin0^^2 0^2 = 008 0^2- (b) Tbe ^-dependence of 05^2 

and ac'i 2 - The asymptotic value of 05^2 is dm|^Si 2 Ci 2 /cf 3 « 0.014 = 0{e‘^\Sm‘^i\). 
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Figure 21. The dependence of (a) (j) and (b) 6*33 = 0 i2,+(j) on /3 = — log^ (a/ldmgil) for the normal 
(NH) and inverted (IH) mass hierarchies. 


where = cos (/>, = sin cf), and 

2dC32Cl3Si3 


tan 2 (j) = — 


d sin 2013 


dmgi — dsfg — A_ 


— 5 m2is\2) + « cos 20i3 


(A.14) 


The angle (j) is in the fourth quadrant when dm^i > 0 (normal hierarchy), and in the 
hrst quadrant when < 0 (inverted hierarchy). The /3-dependence of (j) is shown 
in Fig. 21(a) for both mass hierarchies. 

Using W, we obtain 

77o = W^h'qW 

Azp 01'5 i2C13S13'S(/) 

— aSi2Cl3Sl3'S0 A_,_ 

0 —dSi2Cl3'Sl3C</, 


0 

-dSi2Ci3Si3C,ji 

a: 


(A.15) 
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(a) (b) 

Figure 22. The /3-dependence of Aj_ for the (a) normal and (b) inverted mass hierarchies. 


where the upper signs are for the > 0 case and the lower signs are for the 

< 0 case, with 




[((5m|i - dsfg) -b A_] ± J- asfg) - A_]2 -p 4(oc;2Si3Ci3)^ 


(A.16) 


As /3 is increased beyond 0, the A^ asymptote to 

,2 „2 , ;:_2 „2 „2 
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(3m3iCi3 + (5m2iSi2Si3 , 

-a + + 5m2iS?2C?3 > 


(A.17) 


for both mass hierarchies. Note that A_ < 0 < a!^ for the > 0 case, while 

both A^ < 0 for the Sm^i < 0 case. The /3-dependences of A^ are shown in Fig. 22. 
Order-of-magnitude-wise, we have 
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(A.19) 


For the off-diagonal terms, we have dsg 2 = ci 3 = 0(1), S 13 = 0(e), and 


S 0 = 0 (e), = 0 ( 1 ), if 6 mli > 0 , 
S 0 = 0 ( 1 ) , ^ = 0 (e) , if 6 mli < 0 . 


(A.20) 


•/// . 


Thus, looking at the elements of Hq in that range, we find: 

' 0 {e-^/l) O(eVe^) 0 
H'o = \6ml\ 0(eVe3) 0(e2) 0(eVe") 

0 0(e3/e^) 0(l/e-^) 


(A.21) 
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where the elements with two entries denote the two different mass hierarchies, 0(NH/IH), 
and we see that further diagonalization require angle of order O(e^). Therefore, Hq 
is approximately diagonal. 


A.2.3 7 ^ 0 Case, Third Rotation 

Next, we consider the t] ^ 0 case. Performing the same (1,2) rotation V on h'^ = h'q — 
arjM'^ as we did on h'q, the part is transformed to 


V^M'^{ei2,Bis,923,S)V 
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(A.22) 


Using 9 i2 0, ac'i 2 —> a, asi 2 —^ 0{€‘^)\5rn^i\ as f3 is increased beyond —2, we can 

approximate 
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Performing the (1,3) rotation W next, we find 
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(A.24) 


where s'^g = sin 0 ^ 3 , c'j^g = cos 0 ^ 3 . The angle = 0 i 3 + cj) can be calculated directly 
without the need to calculate cf) using 


tan 20^3 = 


— 517121812) sin ( 20 i 3 ) 


{5m 


31 - ^"^21® 


2 i 

12J 


cos(20i3 ) + a 


(A.25) 


and its /3-dependence is shown in Fig. 21(b). In contrast to the neutrino case, 0^3 increases 
rapidly to 7 r /2 when 5m‘^i < 0, while damping to zero when 5m^i > 0 once /3 is increased 
above zero. Consequently, d cos 0^3 for the 5m‘^i < 0 case, and dsin 0]^3 for the 5m^i > 0 
case plateau to 013513(1 — €^s\2)\5m\i\ = 0{e)\5m\i\ as /3 is increased as shown in Fig. 23. 
Note that in the 5m\i < 0 case, dcos 0 i 3 increased to 0{l)\5m\i\ in the vicinity of /3 = 0 
before plateauing to 0{e)\5m'^^\. As in the neutrino case with 5m‘^i > 0, this will cause 
a slight problem with our approximation later. We now look at the normal and inverted 
mass hierarchy cases separately. 
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Figure 23. The dependence of sin^'^g and cosS^g on /3 = — log^ {a/\5m\P\) for the (a) normal and 
(b) inverted mass hierarchies. The dependence of asin0gg and acos9[^ on /3 = — log^ {a/\5m\P\) 
for the (c) normal and (d) inverted mass hierarchies. 


1. Srrt^^ > 0 Case 

In the > 0 case dsgg —)• 0{e)\6m‘li\ as /3 is increased beyond 0. Therefore, we 

can approximate 
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where we have dropped off-diagonal terms of order 0{e^)\5moi\ or smaller. Define 
the matrix X as 
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Figure 24. Normal hierarchy case. The /3-dependence of x and for several values of rj with 
■S23 = 0-4. 
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Note that 


- ^*23 < X < 0 for 7/ > 0 , 0 < X < ^ — ^23 for 7? < 0 . (A.29) 


The /3-dependence of x is shown in Fig. 24 for several values of r], both positive 
(Fig. 24(a)) and negative (Fig. 24(b)). 


Using X, we find 
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where 
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^ (A.31) 


Thus, is approximately diagonal. The asymptotic forms of Xx± at /3 S> 0 are 
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(A.32) 


Note that have the same asymptotics as Xx± fo'- neutrino case except with 
the 7 / > 0 and rj < 0 cases reversed. The /3-dependence of X'x± are shown in Figs. 24(c) 
to (f). 


2. 5rn^^ < 0 Case 

In the 5777,33 < 0 case 0033 —)• 0(e) 15777 , 33 ! as /3 is increased beyond 0. Therefore, we 
can approximate 
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(A.33) 


where we have dropped off-diagonal terms of order 0 (e^)| 5777 | 3 | or smaller. Define 
the matrix Y as 

C-'ip 


Y = 


0 


-s^e 0 


0 


0 1 


(A.34) 
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Figure 25. Inverted hierarchy case. The /3-dependence of ip and for several values of rj with 
■S23 = 0-4. 


where 


cosip, = sin-i/), and 


tan 27/3 = 


2ar]Si^ sin(2023) 

(A+ - A+) -b 077(1 -b sfg) 008(2023) 

2ar] sin(2023) 

[(Jmlicfs - - 'S?2Si3)] + 207/cos(2023) ' 


(A.35) 
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Note that 


0 < tjj < ^ — 023 for r/ > 0 , —023 < tp < 0 for r/ < 0 . (A.36) 


The /3-dependence of ip is shown in Fig. 25 for several values of rj, both positive 
(Fig. 25(a)) and negative (Fig. 25(b)). 


Using Y, we find 



y^h'''y 


'x'y- 0 0 ‘ 

0 Ay+ 0 


(A.37) 


where 


Ay± 

(A+ -I- A+ - aj/c)! cos 2023) ± \l [A+ - A+ -|- 07/(1 -|- s)|) cos 26 * 23 ]^ + ^(ar/s'^g sin 26 * 23 )^ 

^ 2 ■ 

(A.38) 

Therefore, h'”'_ is approximately diagonal. The asymptotic forms of Ayj_ at /3 2> 0 
are 


Y" -1 I 

Ay_|_ —>■ a|?7| -|- 

Y'" ^ I I 

Ay_ —> + 


-1 1^13^23 + '5w2i(Ci2C23 + S12S13S23) 

-|^”^3i|c? 3«23 + <5"lil(c?2ci3 + S?2«13«23) 

-|(5m|i|c?3ci3 + 


for r] > 0 
for r] < 0 

for ?7 > 0 
for ?7 < 0 


(A.39) 


Note that Ay^ have the same asymptotics as Ay^ for the neutrino case except with 
the 7 / > 0 and i] < 0 cases reversed. The /3-dependence of ip and Ay^ are shown in 
Fig. 25. 


A.3 Effective Mixing Angles for Anti-Neutrinos 

We have discovered that the unitary matrix which approximately diagonalizes H is U* = 
U*QIVW X when > 0, and U* = U*Q^VWY when Sin'll < 0. Taking the complex 
conjugate, these are respectively IJ = UQ 3 VW X* when 6 m^i > 0, and U = UQ 3 VWY* 
when < 0 . 

In the following, we rewrite the mixing matrix in matter U into the form 

U = R23{023,O)Rl3{013,mi2(012,O) , (A.40) 

absorbing the extra mixing angles and CP phase into appropriate definitions of the ‘running’ 
parameters @ 12 , ^ 13 , 023, and 5. As in the neutrino case, frequent use is made of Eq. (3.54). 

• > 0 Case: 

Using Eq. (3.54), it is straightforward to show that 

U = UQsVWX* 
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= -^23(^*235 0)i?i3(6*13, <5)-Rl2(6*12, 0 ) Q 3 i?l 2 (¥^, 0 ) i?l3((/>, 0 ) i?23(X) “<5) 

'-V-' '-V--V--v-" 

U V W X* 

= J^23(^23, O)Q3Ri3(013, O)i?i2(012, 0)i?i2(^, O)iii3(0, 0)i?23(X) “<5) 

= R23{(^23, 0 )QsRi 3 { 6 is, 0 )Ri 2 { 0 i 2 0 )Ri 3 {(j), 0)i?23(X) “^) 

= ^12 

= -^ 23(^*235 0)(53-Ri3(^13) O)-Ri2(012, 0)-Ri 3((/), 0)i?23(X) —<^) ; (A.41) 

where in the last and penultimate lines we have combined the two 12 -rotations into 
one. We now commute i?i 3 ((/i, 0)i?23(X) through the other mixing matrices to the 
left as follows: 


Step 1 : Commutation of iii 3 (()), 0) through Ri2{0i2,0). 

In the range /3 > 0, the angle 0[2 is approximately equal to zero, so we can 
approximate 

■"i 0 0 ’ 

Ruie^.O) ^ Ri2{0,0) = 


0 1 0 
0 0 1 


Note that 


f?12(0, 0)i2i3((/), 0) = i?i3((?!), 0)i?i2(0, 0) 


(A.42) 


(A.43) 


for any cj). On the other hand, in the range — 1 the angle (p is negligibly 
small, so we can approximate 


-Ri 3(</’;0) Ri i?i3(0,0) 


1 0 0 
0 1 0 
0 0 1 


Note that 

.Ri2(^12) 0)-Ri 3(0, 0) = i?i3(0, O)i?i2(0i2,0) 
for any 9^2 ■ Therefore, for all (3 we have 


.Ri2(6*12! 0)-Ri 3(</’) 0) PS i?i3(0, O)i2i2(0i2, 0) , 


(A.44) 


(A.45) 


(A.46) 


and 


U = -R23(^*23, O)Q3i?i3(0i3, O)i?i2(0i2, 0)i?i3(()), 0)i?23(X) —^) 
~ R23{023, 0)Q3Ri3{Gi3, 0)i?l3((/), O)i?i2(0i2,0)i?23(X) “ 

= R23{023, 0)Q3Ri3{Gi3 + 4>, O)i?l2(0i2, 0)i?23(X) —^) 



= R23{023, m3Rl3{e'i^, 0)i2i2(4,0)i?23(X, -<5) , (A.47) 


where in the last and penultimate lines we have combined the two 13-rotations 
into one. The /3-dependence of ^33 was shown in Fig. 21(b). 
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Step 2 : Commutation of i?23(X)<^) through i?i2(0i2)O). 

In the range ^ 0, the angle 0^^2 is approximately equal to zero as we have 
noted above and we have the approximation given in Eq. (A.42). Note that 

i?12(0,0)i223(X)-<J) = -R23 (X)-<^)-Ri2(0,0) (A.48) 


for any x- On the other hand, in the range /3 < 0, the angle x is negligibly small 
so we can approximate 


R23{X,-S) ^ i?23(0,-5) 


1 0 0 
0 1 0 
0 0 1 


Note that 

i?12(4,0)i?23(0,5) = i?23(0,5)i?12(4,0) 
for any 612. Therefore, for all /3 we see that 


^12(4,0)i?23(X,-<5) « i?23(x,-<5)iil2(4,0) , 


(A.49) 


(A.50) 


(A.51) 


and 


U ^ /223(023,O)Q3^13(4,O)i?23(X,--5)i?12(4,O) • (A.52) 


Step 3 : Commutation of i?23(X) ~^) through 0). 

When Srn^^ > 0 we have O'l^ ~ 0 in the range /3 > 1 so we can approximate 


-Ri 3(^*13!0) i?i3(0,0) 


1 0 0 
0 10 
0 0 1 


(A.53) 


Note that 

i?13(0,0)i223(X)-<J) = -R23 (X)-<^)-Ri3(0,0) (A.54) 

for any x- On the other hand, in the range /3 < 0 the angle x was negligibly 
small so that we had Eq. (A.49). Note that 

Rn(0[3,0)Ri2{0,-6) = i?23(0,-<5)i2i3(4,0) , (A.55) 

for any @ 33 . Therefore, for all /? we see that 

Rl3(e'i3,0)R23{%-6) « i?23(x,--5)i?13(4,0) , (A.56) 

and using Eq. (3.54) we obtain 

U « R23{023,0)Q3R23{X, -5)^13(4> 0):Ri2(^12> 0) 

= R23{023, 0)i?23(x, 0)^13(4> ^2(4, 0)Q3 

= R23{023^j^,O)Rl3{0i3,6)Ri2{0i2,O)Q3 

= 4 
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i?23(4> O)i?13(0l3, 5)i2i2(4,0)Q3 


(A.57) 


where in the last and penultimate lines we have combined the two 23-rotations 
into one. The matrix on the far right can be absorbed into the redefinitions 
of Majorana phases and can be dropped. 


Thus, we find that the effective mixing matrix U in the case > 0 can be 

expressed as Eq. (A.40) with the effective mixing angles and effective CP-violating 
phase given approximately by 


12 " 

^4 

= 012 + ^ 

13 - 


= 013 + ^ 

23 " 

" ^23 

= 6*23 + X 

5 ? 

a 5 . 



(A.58) 


• < 0 Case: 

Using Eq. (3.54), we obtain 

U = UQsVWY* 

= 7223(^23) O)-Ri3(013) <5)7?12(012, 0) Q3 0) Rl3{4>: 0) Rui^P, 5) 

U ^ V W Y* 

= R23{S23Y)Q3Rl3{(^13:O)Rl2{0i2:^)^l3{4>iO)Rl2{'P,^) ■ (A.59) 

We now commute Ri3{(p,0)Ri2{'P,S) through the other mixing matrices to the left 
and re-express U as in Eq. (A.40), absorbing the extra mixing angles and CP phase 
into 012, 013, 023, and 5. The first step is the same as the Sm'^i > 0 case, the only 
difference being the /^-dependence of 0i3, which is also shown in Fig. 21(b). 

— Step 2 : Commutation of Ri 2 {'tjj,S) through i2i2(0i2 ,O)- 

In the range f3 >0 the angle 0^^2 is approximately equal to zero as we have noted 
previously, and we have the approximation given in Eq. (A.42). Note that 

-^ 12 ( 0 , 0)i?i2(V', <J) = 72 i2('!/’, <5)7?i2(0, 0) (A.60) 


for any ip. On the other hand, in the range /? < 0 the angle ip is negligibly small 


so that 


i2i2(V’,5) ~ 7?12(0,5) 


1 0 0 
0 1 0 
0 0 1 


(A.61) 


Note that 

i2i2(4,0)i2i2(0,5) = i2i2(0,5)i2i2(4,0) (A.62) 

for any 9 12 - Therefore, for all (5 we see that 


^12(012, 0)i2l2(V', <5) ~ i2l2(V’, 5)^12(012,0) , (A.63) 
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and 


(A.64) 


U « /?23(023,O)Q3i?13(013,O)i?12(V’,<^)^12(0l2,O) . 


Step 3 : Commutation of Ri 2 {'ip, 6 ) through Ri^{ 9 i 2 , 0 ). 
When Srn^i < 0 we have 0^3 | in the range /3 > 1 so that 


-^13(^13)0) ~ -^ 13 ( 2 ’*^ 


TT 


0 0 1 

0 1 0 

-10 0 


Note that 


(A.65) 


^13(^,0)i?l2(V',<^) = ^23(^,-5)i?i3(^,0) (A. 66 ) 


for all ip. On the other hand, in the range /3 < 0 the angle ip was negligibly 
small so that we had the approximation Eq. (3.75). Note that 


i?13(0i3,O)i?i2(O,,5) = i?23 (0,-5)7213(013,0) 
for all 0 ^ 3 . Therefore, for all /3 we see that 

i?i3(4,0)72i2(^,5) ~ 7223(^,-5)7213(4,0), 


(A.67) 


(A. 68 ) 


and using Eq. (3.54) we obtain 

U « 7223(023, 0)g3i?^3(^, -5)^13(013, 0)^12(012, 0) 

= 7223(023 , 0)7223(V^, O)72i3(0'i3, 5)72i2(0i2, 0)Q3 
= 7223(02^^, O)72i3(0'i3, 5)72i2(0'i2,0)Q3 

= ^23 

= 7223(4, 0)72i3(4, -5)72i2(0'i2, 0)^3 , (A.69) 


where in the last and penultimate lines we have combined the two 23-rotations 
into one. The matrix ^3 on the far right can be absorbed into redefinitions of 
the Majorana phases and can be dropped. 


Thus, we find that the effective mixing matrix U for anti-neutrinos in the case Sm^i < 
0 can be expressed as Eq. (A.40) with the effective mixing angles and effective CP- 
violating phase given approximately by 


012 ~ 012 = 012 + ^ , 

013 « 033 = 9i3 + 4> , 

023 « 023 = 023 + V’ , 


(A.70) 
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p 


p 


(a) Normal Hierarchy 


(b) Inverted Hierarchy 


Figure 26. The /3-dependence of 6*23 for the (a) normal and (b) inverted hierarchies for several 
values of 77 with s^g = 0.4. 


A.4 Summary of Anti-Neutrino Case 

To summarize what we have learned, the inclusion of the dr/M^ term in the effective 
Hamiltonian shifts 623 to 6*23 = O 23 -|- X for the > 0 case, and to ^23 = ^23 + V' for the 
(5m|g < 0 case. For both cases, ^23 can be calculated directly without calculating x or tp 
first via the expression 


tan 2023 


[(5m|icf3 - dm|g(cf2 - ■Si2'Si3)] sin 26*23 
[dm|gCg 3 — ( 5 m 2 g(cg 2 — Si 2 '®i 3 )] cos 26*23 + ‘ 2 ar] 


(A.71) 


vr 


Note that as (5 is increased, ^'23 runs toward — if < 0, while it runs toward 0 if 

> 0. The /3-dependence of 623 is shown in Fig. 26. The CP-violating phase 6 
remains unaffected and maintains its vacuum value. 

The running of the effective mass-squared differences are also modified in the range 
/3 > 0. For the Sm^i > 0 case, A 2 and A 3 show extra running, while for the < 0 case, 
it is Ai and A 2 that show extra running, cf. Figs. 24 and 25. 


A.5 Discussion at the Probability Level 

Now, we demonstrate how these lepton-flavor-conserving NSI parameters affect the anti¬ 
neutrino oscillation probability for the various appearance and the disappearance channels. 
In Fig. 27, we compare our approximate oscillation probabilities (blue curves) as a 

function of the neutrino energy against the exact numerical results (red curves) assuming 
r] = 0.1, C = 0 (left panels) and r/ = —0.1,C = 0 (right panels). The upper (lower) panels 
are drawn for 2300 km (8770 km) baseline. Here, in all the panels, we assume 823 = 40°, 
(5 = 0°, and inverted mass hierarchy {dm^i < 0). We also compare our results with the 
approximate expressions of Asano and Minakata [69] (dashed green curves). 

The accuracy of our analytical approximations as compared to the exact numerical 
results for different values of 623 is shown in Fig. 28. All the plots in Fig. 28 have been 
generated assuming (5 = 0° and inverted mass hierarchy < 0). We take the same 
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L=2300km, (ep^=0.10, err=-O.W)/(//=O.IO, f= 0) L=2300 km, (£^^=-0.10, e„=0.10)/(i7=-0.10, f=0) 




L=8770 km, (£^^=0.10, £„=~0.10)/(;7=0.10, f= 0) L=8770 km, (£^^=-0.10, £„=0.10)/(7;=-0.10, f=0) 




E (GeV) E (GeV) 


Figure 27. ^ Ve transition probability as a function of anti-neutrino energy E in GeV for 

2300 km (8770 km) baseline in upper (lower) panels. We compare the analytical expressions of 
this work and Asano-Minakata [69] against the exact numerical result assuming ry = 0.1, C = 0 (left 
panels) and ry = —0.1, C = 0 (right panels). The solid black curves portray the standard three-flavor 
oscillation probabilities in matter without NSI’s. In all the panels, we consider 6*23 = 40°), 8 = 0°, 
and inverted mass hierarchy. 


choices of ry and C like in Fig. 27 and results are given for 2300 km (upper panels) and 8770 
km (lower panels) baselines. We find that our approximation provides satisfactory match 
with exact numerical results for different values of 023- Fig- 29 presents a comparison of 
our approximate probability expressions (solid curves) against the exact numerical results 
(dashed curves) assuming four different values of the CP-violating phase b at 2300 km. 
Here, we take 023 = 40° and brri^^ < 0. 

In Fig. 30, we depict the —?■ survival probability in the presence of NSI for 
two different values of 023 (40° and 50°) at 8770 km baseline. We present the matching 
between the analytical and numerical results assuming ry = 0.1, C = 0 (left panel) and 
ry = —0.1,C = 0 (right panel). In both the panels, we consider <5 = 0° and < 0. 

Fig. 30 portrays that our approximate expressions match quite nicely with the numerical 
results. 
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L=2300 km, (£^^=0.10, e„=-0.10)/(7;=0.10, (= 0) 



E (GeV) 


L=2300 km, (£^^=-0.10, £rr=0.10)/(7/=-0.10, f= 0) 



E (GeV) 


L=8770 km, (£^^=0.10, £„=-0.10)/(7;=0.10, 0) L=8770 km, (£^^=-0.10, £„=0.10)/(;7=-0.10, (= 0) 



Figure 28. Comparison of our analytical expressions (solid curves) to the exact numerical results 
(dashed curves) in case of anti-neutrino for various values of 6 * 23 ■ We assume (5 = 0° and Smh < 0. 
Upper (lower) panels are for 2300 km (8770 km) baseline. 


L=2300 km, (£^^=0.10, £„=-0.10)/(7;=0.10, (= 0) L=2300 km, (£^^=-0.10, £„=0.10)/(J7=-0.10, (= 0) 




E (GeV) E (GeV) 

Figure 29. Comparison of our analytical expressions (solid curves) to the exact numerical results 
(dashed curves) in case of anti-neutrino for different values of the CP-violating phase S at 2300 km. 
To generate these plots, we consider 623 = 40° and Smh < 0. 
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L=8770 km, (6^^=0.10, erp=-0.10)/(7/=0.10, 0) 



L=8770 km, (€^^=-0.10, e^p=0.10)/(?7=-0.10, 0) 



E (GeV) 


(a) rj = 0.1, C = 0 


(b) rj = -0.1, C = 0 


Figure 30. — >■ survival probability as a function of anti-neutrino energy E in GeV for two 

different values of 023 at 8770 km baseline. Comparison between the analytical and numerical results 
assuming r] = 0.1, ^ = 0 (left panel) and rj = — 0.1,C = 0 (right panel). The standard three-flavor 
oscillation probabilities in matter without NSI’s are also shown. In both the panels, we assume 
i5 = 0° and 5m^i < 0. 


B Comparing Probabilities with Constant Varying Earth Density Pro¬ 
file 

So far, we considered the line-averaged constant Earth matter density for a given baseline 
which has been estimated using the PREM profile to present our results. Now, it would 
be quite interesting to study how the exact numerical probabilities would be affected if 
we consider the more realistic varying Earth density profile instead of the line-averaged 
constant matter density for the baselines as large as 8770 km and 10000 km. In Eig. 31, we 
show the exact numerical probabilities considering the constant and varying Earth density 
profiles for 8700 km (upper panels) and 10000 km (lower panels) baselines. In the figure 
legends, the line-averaged constant matter density cases are denoted by ‘pavg’ and the 
varying Earth density cases are labelled by ‘PREM’. We perform these comparisons for 
both the SM and NSI scenarios. For the NSI’s, we take rj = 0.1, C = 0 (left panels), and 
rj = —0.1,C = 0 (right panels). For the sake of completeness, we also plot the approximate 
probability expressions in the presence of NSI that we derived in this paper assuming 
the line-averaged constant Earth matter density based on the PREM profile. Fig. 31 
clearly shows that though the line-averaged constant matter density probability does not 
completely overlap with the PREM-based varying matter density profile probability, it is 
still fairly accurate. Moreover, the effect due to the inclusion of flavor-diagonal NSI’s is 
correctly captured in our approximate analytical expressions. 
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L=8770 km, (£^^=0.10, e„=-0.10)/(;;=0.10, (= 0) 


L=8770 km, (£^^=-0.10, £„=+0.10)/(7;=-0.10, (= 0) 




L=10000 km, (£p^=0.10, £„=-0.10)/(;;=0.10, (= 0) 




2 4 6 8 10 12 


E (GeV) 


Figure 31. ^ i/e transition probability as a function of neutrino energy E in GeV for 8770 

km (10000 km) baseline in upper (lower) panels. We compare the exact numerical probabilities 
considering the constant (denoted by ‘pavg’)i varying (labelled by ‘PREM’) Earth density 
profiles for both the SM and NSI cases. We also plot our approximate analytical expressions in the 
presence of NSL For the NSI’s, we take 77 = 0.1, ^ = 0 (left panels), and 77 = —0.1,C = 0 (right 
panels). In all the panels, we consider 023 = 40°, i5 = 0°, and normal mass hierarchy. 
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